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IMPROVEMENTS  IN  EMPIRICAL  MODELLING 
OP  THE  WORLD-WIDE  IONOSPHERE 


1 .  INTRODUCTION 

Numerical  Models  of  the  ionosphere  are  needed  for  very  diffe¬ 
rent  applications:  radio  wave  propagation,  environaental  stu¬ 
dies,  estimates  for  theoretical  considerations  are  a  few  of 
these.  While  for  experiments  with  radio  waves  it  is  important 
to  have  a  rather  accurate  vertical  electron  density  profile  at 
the  Interesting  position,  predictions  of  radio  wave  propaga¬ 
tion  often  need  this  knowledge  over  a  larger  sector  of  Earth. 
For  almost  all  radio  problems  the  electrons  alone  are  of  impor¬ 
tance,  such  that  the  composition  of  the  heavier  ions  must  not 
be  known;  also,  in  this  context,  the  electron  and  ion  tempera¬ 
tures  are  not  important.  These  parameters  are,  however,  of 
great  interest  for  environaental  studies  and  for  checking 
aeronoaic  theories. 

1.1  The  CCIR  Peak  Program.  Therefore,  according  to  their 
uses,  there  exist  different  models.  As  for  radio  wave  propa¬ 
gation  in  the  hf  range,  the  Coaiti  Consul tatif  International 
dee  Radlocoamunications  (C.C.I.R.)  has  accepted  in  1967  a  com¬ 
puter  program1  specifying  monthly  median  peak  data  of  the  iono¬ 
sphere,  namely:  foF2,  the  critical  frequency  of  the  F2-layer 
and  M(3ooo)F2,  a  propagation  parameter  used  for  computing  the 
maximum  usable  frequency  (MUF)  for  one-hope  ionospheric  propa¬ 
gation  over  a  given  distance.  Both  parameters  are  easily  deri- 


C.C.I.R.,  Atlas  of  ionospheric  characteristics,  Rept.  No.34o, 
Union  Internationale  des  Tllicoaaunications,  Genlve  1983 
(original  1967/74). 
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▼ed  from  (bottom-side)  ionograas  .  Proa  the  square  of  foF2  one 
obtains  directly  the  peak  electron  density,  NaF2.  M(3ooo)F2  is 
narrowly  connected  with  the  peak  height,  haF2.  Different  con¬ 
necting  relations  hare  been  indicated  in  the  literature***^. 

The  CCIR  peak  prograa^  derives  entirely  froa  a  large  set  of 
aonthly  aean  data  aeasured  at  about  loo  ionospheric  sounding 
stations  all  over  the  world  and  for  aany  years.  The  geographic 
distribution  of  the  stations  is  not  uniform,  since  there  were 
auch  less  stations  in  the  southern  heaisphere  than  in  the  nor¬ 
thern  one  and  alaost  none  in  the  oceans.  The  method  of  repre- 


2 

V.  Roy  Piggott  and  Karl  M.A.  Rawer,  U.R.S.I.  Handbook  of 
Ionograa  Interpretation  and  Reduction,  Elsevier,  Amsterdam 
1961  (and  translations  into  french,  Japanese  and  russian 
languages).  Seed,  edition:  Rept.  UAG-23,  World  Data  Center  A 
S.T.P.,  Boulder,  Co.,  U.S.A.  1972.  Revision  of  chapts.  1-4 
as  Rent.  UAG-23 A.  ibidea  1978. 

^  T.  Shiaasaki,  World-wide  daily  variations  in  the  height  of 
the  aaxiaua  electron  density  of  the  ionospheric  F2-layer, 

J.  Radio  Res.  Labs.  (Tokyo)  2,  83-97  (1939). 

4  R.B.  Bent  and  S.K.  Llewellyn,  Description  of  the  1963-71 
ionospheric  model  used  in  the  Definitive  System  (DODS),  Rent. 
DBA- Systems.  Melburae,  FI.,  U.S.A.  197o. 

P.A.  Bradley  and  J.R.  Dudeney,  Vertical  distribution  of  elec¬ 
tron  concentration  in  the  ionosphere,  J.  Ataos.  Terr.  Phys.  ££ 
2131-2146  (1973). 

^  D.  Bilitsa  and  R.  Eyfrig,  Modell  sur  Darstellung  der  Hohe  des 
F2-Maxiauas  alt  Hilfe  des  N(3ooo)F2-Wertes  des  CCIR,  Klein- 
heubacher  Berichte  21.,  167-174  (1978). 

*  C.C.l.R. ,  Atlas  of  ionospheric  characteristics,  Rept.  No.  34o, 
Union  Internationale  des  Tfclicomaunications,  Genive  1983. 
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sentation,  due  to  Jones  and  Gallet^,  starts  with  a  Fourier  ana¬ 
lysis  (of  seventh  order)  of  the  (aonthly  Bean)  diurnal  varia¬ 
tion.  It  ends  up  with  13  Fourier-coefficients  for  each  station. 
Legendre  (spherical  functions)  analysis  is  then  applied  separa¬ 
tely  to  each  of  these  coefficients,  thus  establishing  13 
world-wide  naps.  None  of  these  has  a  direct  geophysical  mea- 
ning,  but  all  together  with  the  relevant  sine  and  cosine  func¬ 
tions  can  be  used  to  coapose  a  world-wide  aap. 

This  rather  coaplicated  procedure,  introduced  after  trials  had 
failed  to  apply  the  Legendre  procedure  directly  to  foF2  or 
M(3ooo)F2.  Such  representation  has  an  understandable  tendency 
to  saooth-out  the  very  characteristic  steep  slopes  which  occur, 
in  particular  after  sunrise. 

The  fact  that  the  world-wide  distribution  of  the  baBic  inputs 
is  far  froa  being  unifora  forced  Gallet  and  Jones  to  develope 
a  special  procedure  for  deteraining  the  Legendre  coefficients. 
(The  traditional  procedure  needs  a  uniform  grid).  However,  when 
they  applied  their  aethod  to  the  aeasured  data,  they  got  nega¬ 
tive  values  of  foF2  in  the  Pacific  (where  their  grid  width  was 
particularly  large).  In  order  to  exclude  such  aberrations, 
they  introduced  so-called  'screen  stations'  by  shifting  coastal 
stations  over  soae  distance  into  the  oceans,  thus  using  the 
saae  data  input  twice.  While  this  trick  aade  the  result  aore 
reasonable,  it  had  a  disappointing  effect  on  the  latitudinal 
variation  over  the  oceans. 

Feeling  that  shifting  along  circles  of  constant  latitude  was 
not  appropriate  because  the  upper  ionosphere  is  under  strong 
geoaagnetlc  control,  so  that  Rawer®  proposed  another  latitudi¬ 
nal  coordinate,  now  called  MODIP  (modified  dip)  which  better 


^W.B.  Jones  and  R.M.  Gallet,  Telecom.  J.  2£,  129  (1962); 
ibidea  £2,  18  (1965). 

g 

Karl  M.A.  Rawer,  F2-layer  ionisation,  in:  B.  Landmark  (ed.), 
Advances  in  Upper  Atmosphere  Research,  Pergamon,  Oxford  1963, 
pp.  159-2o7. 


takes  account  of  this  geophysical  conditions.  The  CCIR  program 
derives  fro®  an  analysis  nade  with  this  coordinate  (instead  of 
earlier  geographic  latitude). 

1.2  The  Profile  Shape.  Apart  from  this  peak  program,  in 
view  of  absorption  computations,  CCIR  accepted  later  a  standard 
vertical  profile  following  a  proposal  by  Bradley  and  Dudeney5,9 
This  is  a  combination  of  two  parabolas,  one  for  the  E-  and  one 
for  the  F-bottomside,  linked  by  a  linear  transition  range.  Ano¬ 
ther  vertical  profile  model  is  used  in  the  Bent-model^  which  is 
applied  in  NASA  practice  for  computing  the  different  effects  of 
ionospheric  refraction  on  Earth-satellite  radio  links  in  the 
vhf  and  uhf  ranges.  These  authors  have  a  fourth  order  parabola 
for  the  bottom-side  of  the  ionosphere,  a  second  order  parabola 
for  the  topside  which  on  top  is  linked  successively  with  three 
exponential  sections.  Since  the  bottomside  enters  into  the  com¬ 
putations  only  with  its  total  electron  content,  the  lower  iono¬ 
sphere  is  sumaarisingly  taken  account  of  by  an  increased  thick¬ 
ness  of  the  fourth  order  bottomside  parabola.  The  Bent  model  is 
particularly  important  by  the  fact  that  it  is  the  only  readily 
accessible  summary  of  a  very  large  set  of  about  lo.ooo  topside 
profiles  derived  from  ionograms  gathered  by  the  ALOUETTE  satel¬ 
lites.  All  these  records  are  from  receiving  stations  at  diffe¬ 
rent  latitudes,  but  all  in  the  longitude  range  of  the  American 
continent.  The  latitudinal  variation  in  the  Bent-model  is  dis¬ 
continuous,  admitting  only  three  latitude  ranges.  It  is  unfor¬ 
tunate  that  the  whole  equatorial  range  (+  3o°  geomagnetic)  was 
thrown  into  just  one  latitude  range.  Thus,  in  the  equatorial 
cone  Bent's  reproduction  of  the  outer  ionosphere  profile  can¬ 
not  show  the  important  variations  with  latitude  which  occur 
Inside  this  zone. 

5  P.A.  Bradley  and  J.R.  Dudeney,  Vertical  distribution  of  elec¬ 
tron  concentration  in  the  ionosphere,  J.  Atmos.  Terr.  Phys. 
2131-2146  (1973). 

9  J.R.  Dudeney,  J.  Atmos.  Terr.  Phys.  40,  195  (1978). 

*  R.B.  Bent  and  S.K.  Llewellyn,  Description  of  the  1965-1971 
ionospheric  model  used  in  the  Definitive  System  (DODS),  Rept, 
DBA- Systems ,  Melbourne,  PI.,  U.S.A.  197o. 


1 .3  The  International  Reference  Ionosphere. 

The  International  Reference  Ionosphere  (IRI)  is  an  internatio¬ 
nal  project  which  originated  in  COSPAR  (the  Space  Research  Com- 
■ittee  of  ICSU)  and  is  co-sponsored  by  the  International  Union 
of  Radio  Science,  URSI.  In  view  of  enTironaental  applications, 
which  are  of  first  interest  to  COSPAR,  apart  from  electron  den¬ 
sities  this  aodel  gives  also  electron  and  ion  teaperatures  and 
the  cheaical  coaposition  of  the  positive  ion  population.  The 
first  IRI  was  published  in  197810;  a  later  publication  called 
IRI-7911  corresponds,  in  fact,  to  the  state  of  198o.  The  aodel 
is  fully  coaputerised.  During  the  years,  the  'URSI-COSPAR  Task 
Group  on  the  IRI*  has  introduced  different  improvements  and  the 
coaputer  program  was,  accordingly,  refurbished  several  times. 
The  actually  valid  version  is  IRI-9. 

As  for  the  main  peak,  the  IRI  electron  density  profile,  Ne(s), 
depends  entirely  on  the  CCIR  program1.  So,  the  profile  is  given 
in  reduced  fora  coaparing  with  the  ?2  peak  density,  NmF2 .  The 
logaritha  of  Ne/NaP2  is  described  layerwise  by  aatheaatical  ex¬ 
pressions  which  are  smoothly  connected  at  the  interfaces.  This 
is  a  rather  involved  procedure,  but  it  has  the  advantage  that 
measured  characteristics  like  foE,  foPI  etc.  can  easily  be  in¬ 
troduced. 

10  K.  Rawer,  S.  Ramakrishnan  and  D.  Eilitsa,  International  Re¬ 
ference  Ionosphere  1978,  International  Union  of  Radio  Scien¬ 
ce  (URSI),  Brussels,  B.,  1978. 

11  K.  Rawer  (chan.),  J.V.  Lincoln  and  R.O.  ConKright  (eds.). 
International  Reference  Ionosphere  -  IRI  79 ,  World  Data 
Center  A  (S.T.P.),  Boulder,  Co.,  U.S.A.  1981. 

1  C.C.I.R.,  Atlas  of  ionospheric  characteristics,  Rept.  No. 

21®,  Union  Internationale  des  Telecommunications,  Geneve 
1983. 
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1.4  Booker *8  Proposal.  Right  at  the  beginning  of  the  IRI 
studies,  Booker12  proposed  another  reproduction  of  the  He- 
profile  which  is  fully  analytic.  It  is  apparent  that  such 
representation  is  advantageous  for  ware  propagation  conputa- 
tions.  Booker  starts  froa  a  graph  of  the  logarithaic  deriva¬ 
tive  of  the  profile,  d  log  Ne/dz,  which  he  approxiaates  step¬ 
wise  by  constant  values.  In  order  to  aake  this  'skeleton  func¬ 
tion'  analytic,  Booker  replaces  each  of  the  discontinuous  steps 
by  one  fully  continuous  Bpsteln-step-function.  Eps0,  which  is 
given  by 

Bps0  (z;HX,SC)  -  1/(1  ♦  exp(-  |))  (1) 

with  (z-HX)/SC. 

It  is  evident  that,  aeyaptotically,  EpsQ  approaches  one  at  the 
right  hand  side  and  zero  at  left.  The  aain  variation  occurs  in¬ 
side  the  width  SC  and  is  centered  at  HX,  where  the  function  ta¬ 
kes  the  value  1/2.  By  linearly  coabining  several  functions  Eps0 
with  different,  suitably  chosen  paraaeters  (plus  a  constant), 
the  original  skeleton  can  quite  well  be  approxiaated .  The  pro¬ 
file  itself  aust  then  be  found  by  integration.  The  integral 

13 

function  to  Eps0  reads 

Eps.l  (z;HX,SC)  »  ln(1  «■  exp(  J  ))  (2) 

and  is  called  Epstein-transition .  It  has  for  right  hand  asynp- 
tote  the  linear  function  y  ■  f  with  slope  1 ,  while  the  left 
hand  asyaptote  is  the  e-axis  (as  for  £ps0).  For  coapleteness, 
we  yet  noie  the  derivative  of  EpsQ  which  reads 


12 

H.6.  Booker,  Fitting  of  aulti-region  ionospheric  profiles  of 
electron  density  by  a  single  analytic  function  of  height, 

J.  Ataos.  Terr.  Phys.  22'  619-623  (1977). 

^  K.  Rawer,  Replaceaent  of  the  present  sub-peak  plasma  density 
profile  by  a  unique  expression,  Adv.  Space  Res.  2,  No.  lo, 
183-19o  (1982).  (Printing  error  in  formula  for  £ps^). 
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Epe^BjHX.SC)  -  exp(/)/(1  ♦  exp(f  ))2.  (3) 

It  peaks  at  0,  i.e.  s  ■  HI,  and  has  on  both  sides  the  s- 
axis  as  asymptote.  It  is,  therfore,  called  Epstein-layer. 


After  Integration,  Booker's  reproduction  of  the  profile  is  a 
sua  of  Eps_.,  functions  plus  a  linear  tern  (which  etens  fron  the 
constant  in  the  skeleton).  His  fitting  procedure  is  essentially 
layer-by-layer,  i.e.  by  successive  approximation. 

Booker' s  proposal  was  discussed  in  the  Task  Group  on  IRI.  Since 
in  practice  one  has  sons  characteristic  values  rather  than  the 
full  derivative  profile,  it  appeared  difficult  to  meet  these 
constraints  when  following  Booker's  proposal.  It  appeared  at 
that  time  that  a  lengthy  trial-and-error  process  applied  to 
each  Individual  profile  was  the  only  way  to  apply  the  Booker 
method  with  given  constraints.  A  computer  program  going  that 
way  was  besides  presented  in  the  first  (1978)  edition  of  IRI, 
but  was  not  often  used. 

1 .5  The  Topside  Profile.  In  the  same  edition,  and  in  all 
that  followed,  a  description  in  terms  of  two  Epstein  transition 
functions  (plus  a  linear  one)  was,  however,  applied  to  the  top¬ 
side.  This  description  (which  is  due  to  S.  Ramakrishnan )  was 
indirectly  based  on  the  data  set  gathered  by  the  ALOUETTE  sa¬ 
tellites  since  the  descriptive  function  was  fitted  to  Bent's 
description.  Although  this  latter  is  discontinuous,  the  new 
description  was  made  continuous  in  all  parameters  (latitude, 
foF2,  solar  activity). 

The  IRI  topside  formula  depends  on  eight  parameters,  namely  one 
HX,  one  SC  and  one  amplitude  to  each  of  the  two  Eps_.|  members 
plus  the  two  parameters  of  the  additional  linear  function. 

These  were  fitted  so  as  to  have  a  peak  at  the  given  altitude 
hmF2.  With  the  first  Eps_1t  i.e.  with  HX1,  SCI  and  A1  Ramakrish¬ 
nan  approached  the  given  Bent-profile  in  the  loo  to  2oo  km 


height  range,  just  above  the  peak  and  with  HX2,  SC2  and  A2,  he 
adapted  the  function  to  the  Bent-profile  near  7oo  to  800  km  of 
altitude.  This  was  in  agreement  with  Booker's  proposal.  How¬ 
ever,  in  order  to  improve  the  description,  a  final  correction 
was  made  taking  simultaneously  account  of  both  functions  to¬ 
gether. 

It  was  evident  that  a  generalisation  of  this  latter  procedure 
to  the  whole  ionospheric  profile  was  not  very  promising,  be¬ 
cause  one  had  then  to  deal  with  a  considerable  number  of  Eps_., 
members  (between  seven  and  twelve)  and,  consequently,  with 
quite  a  large  number  of  parameters  to  be  determined  at  one  time 
(certainly  more  than  twenty). 

1.6  A  Way  to  Ease  Fitting.  Later,  Rawer  made  another  pro¬ 
posal1*  in  order  to  reduce  the  number  of  parameters  to  be  fitted 
at  one  time.  He  felt  that  by  subdividing  the  height  range  into, 
aay  three  ranges,  one  might  be  able  to  come  to  a  reasonable 
number  of  unknowns  to  be  determined  Independently  in  each  range. 
He  applies  filter  functions  by  which  a  given  height  range  night 
be  selected.  A  function  of  the  following  kind  is  chosen: 

Pi(*;sk_1 ,sk)  -  Eps0(s;»k_1 ,s)  -  Eps0(z;zk,s)  (4) 

where  *k  is  the  upper,  zk_1  the  lower  limit  of  the  selected 
range  and  s  is  a  scale  which  night  be  of  the  order  of  a  few  km. 
Apart  from  small  ranges  near  z*zk-1  and  z-zk,  ?ik  18  e88enti“ 
ally  one  in  the  range  zk-1  <  t  <  zk,  and  zero  outside.  How,  if 
for  this  range  any  description  of  the  profile  has  been  establi¬ 
shed  as  an  anlytic  function,  Ne(z),  this  latter  might  take  va¬ 
lues  far  away  from  the  given  profile  outside  this  range.  How¬ 
ever,  when  multiplying  Ne(z)  by  the  relevant  filter  function 
Fi,  the  product  is  essentially  zero  outside  that  range  and 

**  K.  Rawer,  Analytical  description  of  profiles  through  plane¬ 
tary  atmospheres,  Acta  Astronaut.  JJ_,  6o7-6o8,  1984. 


identical  with  Ne(z)  inside  it.  Applying  this  procedure  sepa¬ 
rately  in  each  range  and  adding  up  the  different  products,  one 
obtains  a  correct  reproduction  of  the  whole  electron  density 
profile,  though  fitting  was  applied  individually  and  indepen¬ 
dently  in  each  range. 

At  the  range  limits  one  has  a  continuous  change-over  from  one 
to  the  next  individual  description.  In  order  to  have  continu¬ 
ity  in  value  and  slope,  it  is  requested  that  at  the  cutting 
height  the  two  individual  functions  take  the  same  value  and 
derivative  from  both  sides.  It  is  proposed  to  have  the  peaks 
of  F2-  and  E-layer  as  range  limits,  such  that  one  has  three 
subranges:  topside  -  middle  ionosphere  -  lower  ionosphere. 

The  values  at  the  limits  are  then  Nm?2  and  NmE,  with  zero 
derivative  at  both  limits. 

1.7  Two  Different  Pitting  Methods.  Booker's  method  of  suc¬ 
cessive  fitting  ('layer-by-layer')  is  rather  easy  to  apply.  It 
has,  however,  an  important  disadvantage,  because  each  £ps_1 
function  describes  the  profile  only  in  a  subrange  within  which 
the  logarithmic  derivative  should  essentially  be  constant.  This 
is  so  because  the  skeleton  function  is  piecewise  approached;  an 
overlap  of  the  individual  Eps0  functions  is  not  allowed  in 
Booker's  method.  It  is  evident  that  then  the  number  of  subran¬ 
ges  and  that  of  individual  functions  might  become  rather  large. 
This  might  be  acceptable  at  profile  analysis  where  it  is  inten¬ 
ded  to  reproduce  just  one  given  profile,  but  it  is  not  suitable 
at  profile  synthesis,  l.e.  for  use  in  a  model. 

The  number  can  be  considerably  reduced  when  overlap  of  the  dif¬ 
ferent  functions  is  admitted.  In  that  case,  however,  the  fit¬ 
ting  can  no  more  be  made  successively,  function  by  function, 
but  has  to  be  made  simultaneously  for  all  functions  which  are 
needed  in  a  given  range.  This  method1^  will  be  applied  in  the 

ie 

K.  Rawer,  L.  Bilitea  and  T.L.  Gulyaeva,  New  formulas  for  the 

IRI  electron  density  profile  in  the  topside  and  middle  iono¬ 
sphere,  Adv.  Space  Res.  £,  No.  7,  3-12  (1963). 


following.  It  was,  in  fact,  applied  for  the  first  tine  with  the 
final  adjustaent  of  the  IRI  topside  foraula. 

1 .8  Siting  the  Main  Peak.  With  our  choice  of  the  range  li- 
aits  at  ha?2  and  ha£,  there  exists  for  each  range  one  peak  at 
one  of  the  range  liaits.  Its  height  and  value  are  a  priori  gi¬ 
ven  by  an  external  input,  e.g.  the  CCIR  prograa  or  a  foraula 
of  geophysical  character.  It  is  iaportant  that  these  values 
are  correctly  reproduced. 

The  value  problea  is  easily  resolved  because  in  each  range  we 
reproduce  a  reduced  (logarithaic)  profile,  naaely  log(Ne(z)/Nem). 

Consequently,  the  profile  function  should  be  zero  at  ha. 
Further,  the  derivative  d  logNe/dz  should  also  be  sero  there. 
In  order  to  fulfill  these  two  conditions,  we  need  two  paraae- 
ters  which  can  be  freely  adapted.  The  linear  aeaber  is  adequate 
to  this  end.  So,  the  peak  condition  aight  be  satisfied  by  sui¬ 
tably  choosing  the  linear  aeaber.  Apparently,  this  can  only  be 
done  after  deteraining  the  paraaeters  of  the  different  Spas1 
functions  because  these  appear  in  the  detailed  foraulation  of 
the  two  above  conditions.  On  the  other  hand,  the  linear  aeaber 
is  of  non-negligible  iaportance  in  the  fitting  procedure  it¬ 
self.  Thus,  an  iterative  coaputation  schedule  would  be  needed. 
Rawer1 *  has  shown  how  such  cuabersoae  procedure  can  be  avoided. 
Instead  of  fulfilling  the  peak  conditions  at  the  end,  this  can 
be  done  before  the  fitting  procedure,  provided  an  individual 
linear  tera  is  added  to  each  aeaber  £psa1  in  such  way  that  each 
so  foraed  function  fulfills  the  peak  conditions  individually. 
This  function  is  called  LAY  and  reads 

LAY(z;HX,SC)  -  Eps_1 (z;HX,SC)  -  Eps^ (HMjHX,SC)  - 

-  (z-HM)  •  Eps0(HM;HX,SC)  .  (5) 

^  X.  Rawer,  Analytical  description  of  profiles  through  plane¬ 
tary  ataospheres,  Acta  Astronaut.  1^,  6o7-6o8,  1984. 
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HM  is  the  given  peak  altitude.  By  definition,  LAY  and  its  deri¬ 
vative  as  well  are  sero  at  z*HM,  not  withstanding  the  values  of 
the  geoaetric  parameters  HX  and  SC.  Figure  1  shows  the  depen¬ 
dence  of  LAY-functions  on  different  HX  and  SC  values.  It  ap¬ 
pears  that  the  aain  curvature  occurs,  as  it  must  be,  near  z»HX. 
If  HX  is  auch  lower  than  HM,  the  function  takes  quite  small  va¬ 
lues  in  between,  though  the  peak  (of  value  zero)  occurs  exactly 
at  z-HM. 

When  analyzing  a  profile  in  terns  of  a  number  NFU  of  LAY-func¬ 
tions,  the  number  of  free  parameters  is  3  • NFU.  Of  these,  a  set 
of  NFU  amplitudes  enter  the  analysis  as  linear  parameters, 
while  twice  that  number  of  geometric  parameters,  HX  and  SC,  go 
in  non-linear ly. 

1 .9  Objective  of  the  Study.  When  ionospheric  electron  densi¬ 
ty  profiles  are  to  be  reproduced  along  the  lines  indicated  a- 
bove,  it  is  important  that  the  geoaetric  parameters,  HX  and  SC, 
of  the  different  LAY-functions  can  be  predicted.  In  order  to 
achieve  this  goal,  representative  values  must  be  derived  from 
a  suitable  data  base.  Since,  in  the  functional  representation 

log(Ne(z)/Hem)  -  ^  A^  •  LAY(z;HXi,SCi)  (6) 

i-1 


the  amplitudes  A,  enter  linearly,  these  can  easily  be  adapted 

A  1  £ 

to  given  constraints  in  a  synthesis  program  .  Not  so,  however, 
the  geometric  parameters  HX^  and  SC^.  These  should  be  known  a 
priori  before  synthesis  is  undertaken,  i.e.  when  the  model  is 
used.  So,  the  emphasis  of  our  study  shall  be  on  these  geoaetric 
paraaeters.  It  is  also  iaportant  that  the  miniaua  number  I  of 
aembers  in  Equation  (6)  (needed  to  obtain  a  usable  representa¬ 
tion)  should  be  discussed  for  the  different  geophysical  situa¬ 
tions  which  are  of  interest  in  our  context. 


^  K.  Rawer,  Determining  electron  density  profiles  for  the  mid¬ 
dle  ionosphere,  Adv.  Space  Res.  No.  1o,  43-49  (1983). 
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2.  DATA  BASE 

2.1  Middle  Ionosphere.  Our  aain  data  base  is  a  large  set  of 
profiles  which  were  computed  from  good  quality  ionograms  recor¬ 
ded  at  Lindau/Harz  (51.62°N,  1o.o9°E)  by  V.  Becker1.  The  final 
profile  is  tabulated  with  high  resolution  (plasaa  frequency 
steps  of  1/1 oo  of  foF2).  Of  course,  the  lowest  part  of  the  pro¬ 
file  is  extrapolated  (by  a  cosine  function)  since  there  were  no 
measured  data  below  a  sounding  frequency  of  1  KHz.  Considering 
this  fact  in  our  evaluation,  we  avoided  the  use  of  height  va¬ 
lues  relative  to  plasaa  frequencies  below  o.5  MHz. 

2 

Becker's  coaputation  method  applies  the  often  used  parabolic 
approximation  for  the  part  of  the  profile  immediately  below  the 
peak.  The  critical  frequency  is  obtained  very  accurately  by  this 
fit  which  is  repeatedly  applied  whenever  needed.  In  the  middle 
part  of  the  profile  he  applies  splining  (generally  of  Becond 
order).  The  reading  frequencies  are  individually  chosen  so  that 
to  each  step  in  frequency  corresponds  roughly  the  same  height 
step.  The  "underlying  ionization"  is  taken  account  of  by  an  as¬ 
sumed  cosine  extrapolation  towards  zero  plasma  frequency. 

Though  Becker  has  a  particular  sub- routine  for  taking  account 
of  an  E-F-valley  (i.e.  a  minimum  of  ionization  above  the  E-peak) 
this  was  not  applied  in  almost  all  of  the  profiles  we  have  re¬ 
ceived.  So,  his  results  are  "lowest  acceptable"  profiles.  This 
means:  in  case  a  valley  was  present,  the  real  heights  are  higher 
than  given  by  his  profile.  The  difference  is  the  valley  width 
Just  above  the  E-layer  critical  frequency,  foE;  it  decreases 


We  are  grateful  to  Dr.  Valter  Becker  for  giving  us  two  magne¬ 
tic  tapes  containing  a  large  number  of  such  profiles. 

V.  Becker,  Die  Bestimmung  der  Feins truktur  der  Ionosphhre  aus 
lonogrammen,  Kleinheubacher  Berichte  37-44  (1969)* 
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with  increasing  frequency.  Nevertheless,  the  error  at  the  F2- 
peak,  i.e.  in  ha?2,  sight  be  about  1o  km. 

Otherwise,  aost  of  the  profiles  on  the  tapes  are  reliable. 

There  is  probably  soae  saoothing  introduced  by  the  second  or¬ 
der  splining,  but  this  is  on  line  with  our  intentions:  most 
original  ionograas  show  saall  deforaations  of  the  trace  which 
are  aostly  due  to  gravity  waves  in  the  upper  ataosphere.  For 
our  purpose,  these  aust  be  considered  as  unwanted  perturbations 
of  the  intended  undisturbed  profile.  So,  the  kind  of  saoothing 
introduced  by  the  splitting  is  just  waht  we  need  in  order  to  get 
representative  profiles . 

In  a  very  large  data  set  as  Becker's  it  is  extremely  difficult 
to  exclude  that  a  saall  nuaber  of  the  results  is  in  error.  We 
have  found  a  few  profiles  with  apparent  errors  in  height  sca¬ 
ling  (e.g.  E-peak  at  5o  or  l5o  ka).  Soae  day-time  profiles  show 
no  E-region  -  probably  because  the  corresponding  trace  was  not 
visible  on  the  ionogram.  Profiles  with  apparent  errors  of  this 
kind  were  omitted  from  our  evaluation. 

The  data  set  covers  day  and  night  about  equally  well.  For  most 
evaluated  dates,  ionograms  had  been  taken  at  hourly  intervals 
or  even  mojpe  often  so  that  the  diurnal  variation  is  covered. 
While  the  density  of  observations  is  high  in  this  respect,  the 
periods  treated  are  not  uniformly  distributed  over  seasons  and 
years.  Two  periods  are  particularly  well  covered,  namely: 

-  summer  (27.6.  -  1.7.)  of  1954 

-  spring  (21.3.  -28.4.)  of  1958. 

Quite  clearly,  1958  was  a  year  of  extremely  high  solar  activity 
while  1954  was  about  minimum.  This  is  not  too  bad  for  our  pre¬ 
sent  purpose,  because  we  can  expect  that  the  range  of  parame¬ 
ters  can  be  estimated  when  comparing  the  two  periods. 

The  total  nuaber  of  points  on  tape  determining  the  profile  is 
too  large  for  our  purpose  of  fitting  the  profile  with  a  small 
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number  of  functions.  Also,  thanks  to  saoo thing,  the  information 
about  the  profile  does  not  correspond  to  such  large  input.  We 
have  aade  the  first  fitting  trials  with  every  second  point, 
l.e.  5o  in  total.  Since  the  coaputer  tiae  increases  considera¬ 
bly  with  the  nuaber  of  points  to  be  fitted,  we  finally  used 
every  fourth  point,  i.e.  25  in  total.  This  nuaber  further  is 
reduced  because  the  lower  frequency  data  (sounding  frequency 
below  o.5  MHs)  are  oaitted.  However,  this  was  done  by  putting 
the  weight  to  aero  so  that  the  nuaber  of  points  in  the  algo- 
ritha  reaained  25  always. 

Our  ala  with  these  data  is  getting  characteristics  for  the  aid- 
die  ionosphere.  So,  sounding  results  referring  to  regions  below 
the  E-peak  are  of  no  interest.  To  take  account  of  this,  the 
weight  individually  attributed  to  a  profile  point  was  cut  below 
1o5  ka  by  an  Epstein-step  filter  centered  at  that  altitude 
(with  scale  s«4  ka).  Since  bottoa-side  ionograas  give  no  data 
froa  above  the  72-peak,  filtering  at  the  upper  end  of  the 
height  scale  was  not  needed. 

2.2  lower  Ionosphere.  The  electron  density  profile  shape  be¬ 
low  the  E-peak  is  rather  siallar  to  that  in  the  aiddle  iono¬ 
sphere  with  the  difference  that  the  profile  is  alaost  always 
aonotonous.  The  occurrence  of  a  ainiaua  seeas  to  be  very  rare, 
but  a  turning  point  around  8o  ka  is  typical.  So,  the  saae  ae- 
thods  of  representation  can  be  applied  as  for  the  aiddle  iono¬ 
sphere  . 


a.  McNamara' s  data  base.  By  very  detailed  literature 
survey,  McNaaara^  was  able  to  produce  a  large  collection  of  low¬ 
er  ionosphere  electron  density  profiles  observed  with  different 


L.P.  McNamara,  Ionospheric  D-region  profile  base,  a  collection 
of  computer  accessible  experimental  profiles  of  the  D-  and 
lower  E-regions,  Rept.  UA6-67,  WDC  A  (S.T.P.),  Boulder,  Co., 
U.S.A.  1978.  /We  thank  for  granting  us  a  magnetic  tape  con¬ 
taining  these  data/. 


aethods.  His  profile  set  contains,  of  course,  data  of  very  dif¬ 
ferent  quality.  The  elder  in-situ  aeasureaents  are  very  doubt¬ 
ful  and  even  the  aore  recent  ones  hare  yet  probleas  with  abso¬ 
lute  calibrations.  This  latter  is,  unfortunately,  often  height- 
dependent  so  that  even  the  relative  profile  shape  is  not  gua- 
ranted . 

Y .?.  Raaanaaurty,  who  has  taken  over  to  check  this  data  set, 
established  several  groups  according  to  the  standard  of  the  ob¬ 
servations.  Since  our  fitting  procedure  depends  strongly  on  the 
peak  value  -  the  E-peak  here  -  a  large  nuaber  of  observations 
is  not  useful,  because  the  E-peak  is  not  covered  and  the  abso¬ 
lute  accuracy  of  the  pretended  density  values  is  very  poor. 
Therefore,  only  a  saall  selection  of  profiles  of  this  collec¬ 
tion  could  be  used  for  our  purpose.  It  was,  of  course,  obtained 
at  very  different  locations  and  tiaes.  There  is,  however,  a 
clear  preference  for  day-tiae  and  northern  aiddle  latitudes. 

Night  observations  are  quite  rftre. 

b.  Thuaba  probe  data.  Subbaraya  et  al.*  have  publi¬ 
shed  a  set  of  rocket  aeasureaents  aade  at  Thuaba  in  southern 
India  (8.53°N),  near  the  aagnetic  equator  (dip  o.4o°S).  The 
data  set  is  auch  aore  hoaogeneous  than  McNaaara's  though,  of 
course,  rather  restricted.  The  aeasureaents  were  all  aade  with 
coaparable  equipaent  and  cover  the  altitude  range  froa  60  (day) 
or  83  (night  and  twilight)  up  to  16o  ka.  For  day-tiae  (near 
noon),  a  total  of  twelve  rocket  ascents  was  available.  Night 
aeasureaents  were  aade  at  six  occasions  and  twilight  ones  at 
seven,  three  in  the  aoraing,  four  in  the  evening.  The  authors 
give  a  aean  electron  density  profile  for  day-tiae.  For  the 
other  three  groups,  we  produced  aedians  froa  the  given  profiles. 

*  B.H.  Subbaraya,  Satya  Prakash  and  S.P.  Gupta,  Electron  densi¬ 
ties  in  the  equatorial  lower  ionosphere  froa  the  Langauir  pro¬ 
be  experiaents  conducted  at  Thuaba  during  the  years  1966-1978, 
Rept.  ISRO-PRL-SR- 15-83,  Indian  Space  Research  Organisation, 
Bangalore,  1983. 
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So,  for  each  of  the  four  groups  we  hare  obtained  a  aore  or 
leas  representative  profile.  The  definition  of  this  latter  is 
best  for  day  and  not  too  bad  for  night.  Understandably,  the 
original  twilight  profiles  show  rather  important  differences 
one  against  the  other.  This  aust  be  accepted,  noting  the  quick 
changes  occurring  during  these  periods.  Nevertheless,  we  feel 
that  these  four  profiles  aight,  provisionally,  be  taken  as  re¬ 
presentative  low-latitude  profiles. 

Since  our  analysis  has  to  assuae  an  E-peak,  we  identified  this 
point  individually  on  the  given  profile  knowing,  of  course, 
the  standard  height  range  of  region  E.  Ey  day,  the  experiaen- 
tal  profiles  between  loo  and  Ho  ka  are  soaetiaes  aonotonous, 
soaetiaes  they  show  shallow  ainiaa,  all  with  rather  saall  vari¬ 
ation  in  this  height  range.  The  nean  day- tine  profile  is  aono¬ 
tonous.  In  individual  night  profiles,  however,  above  a  peak  at 
a  height  between  1o4  and  114  ka  up  to  at  least  155  ka,  there  is 
a  large  range  of  definitively  lower  electron  density.  This  is 
also  shown  in  the  nedian. 

c.  Singer’s  profiles.  In  his  thesis  Singer^  produced 
profiles  which,  by  a  kind  of  trial  and  error  aethods,  were  fit¬ 
ted  with  a  rather  large  nuaber  of  propagation  data.  The  question 
whether  such  'conversion'  is  unanbiguous  has  often  been  discus¬ 
sed.  The  author^  feels  that  his  inputs  are  so  nuaerous  that  the 
final  profile  should  be  quite  well  defined.  While  aedians  of 
ln-situ  observations  necessarily  aust  be  saoother  than  the  in¬ 
dividual  profiles,  Singer's  technique  probably  has  a  tendency 
to  produce  rather  sharp  structures  in  the  profile.  This  seeas 
to  be  particularly  so  for  night  tiae.  The  relevant  Singer  pro¬ 
file  is,  in  fact,  a  profile  with  rather  sharp  steps.  We  cannot 
decide  here,  whether  this  is  a  well-based  result  but  have  soae 
doubts  about  this  particular  profile. 


W.  Singer,  Doktor-Dissertation,  Huaboldt-Universitttt  Berlin, 
1976.  /We  thenk  for  forwarding  us  the  profile  tables/. 


3 


MATHEMATICAL  TOOLS 


We  intend  to  reproduce  a  given  profile,  i.e.  a  set  of  between 
15  and  3o  profile  points  by  a  linear  combination  of  three  or, 
at  most,  four  LAY-functions  /Equation  (6)/.  Each  function  con¬ 
tains  two  geometric  parameters,  HX  and  SC,  and  is  multiplied 
in  the  combination  by  a  coefficient  which  might  be  called  am¬ 
plitude  A.  The  task  is  to  find  an  optimum  set  of  the  (unknown) 
parameters  and  amplitudes. 

3.1  Reduced  Error  Sum.  First  of  all,  we  must  define  a  cri¬ 
terion  for  "best  fitting".  In  the  almost  always  used  method  of 
'least  squares'  one  takes  the  sum  of  the  squared  individual  er¬ 
rors  as  the  quantity  which  should  be  minimized  in  the  procedu¬ 
re.  However,  since  we  attribute  different  weights,  W,  to  the 
points,  these  weights  must  be  used  at  the  summation  of  the 
squared  errors.  So,  we  should  consider 

M  2 

Y2  w<k>  *  (Y<k>  -  »<k» 

k-1 


where  Y  is  the  given  value  (at  point  k)  and  F  is  the  corres¬ 
ponding  value  of  the  fitting  function.  This  expression  is  well 
suited  as  criterion  for  fitting. 

However,  since  the  weight  function  W  gives  little  weight  to 
out-of-range  points  (below  a  lower  limit,  HW),  the  number  of 
the  mainly  contributing  points  is  different  for  different  pro¬ 
files.  For  example,  in  Becker's  data  set  HW«1o5  and  the  reading 
points  are  at  equal  steps  of  the  ratio  f/foF2  where  f  is  plas 
ma  frequency  and  foF2  the  peak  plasma  frequency  of  the  F-region 
Thus,  with  a  small  foF2  a  greater  number  of  points  falls  below 
the  E-peak  than  with  a  large  foF2.  For  this  reason,  we  have  to 
expect  a  smaller  sum  for  small  foF2,  i.e.  low  peak  density. 
Since  we  intend  to  have  a  reliable  measure  giving  comparable 
results  under  different  conditions,  a  suitable  reduction  of  the 
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above  expression  is  needed.  It  can  be  obtained  by  dividing  the 
above  sun  by  the  sun  of  the  Individual  weights  (which  in  the 
Becker  schedule  is  smaller  for  small  foF2).  So,  we  define: 


as  reduced  error  sum.  It  gives  the  average  of  the  weighted 
squared  error  per  point. 

3.2  Problems  with  Optimisation. 

a.  Linear  or  non-linear.  Optimisation  is  quite  easy 
if  the  unknown  parameters  appear  linearly  in  the  fitting  func¬ 
tion  F.  In  that  particular  case,  the  method  of  least  squares 
leads  to  a  system  of  linear  equations  with  the  different  para¬ 
meters  as  unknowns.  Since  the  error  sum  is  differentiated  se¬ 
parately  after  each  parameter,  the  number  of  equations  equals 
the  number  of  unknowns.  The  resolution  is  then  obtained  by 
standard  methods  provided  the  determinant  of  the  coefficient 
matrix  is  different  from  zero.  This  condition  is,  however,  on¬ 
ly  satisfied  when  all  the  different  parameters  are  mutually 
independent.  If  there  exists  a  (linear)  connection  between  two 
of  them,  the  determinant  becomes  zero  such  that  no  unique  solu¬ 
tion  can  be  found. 

In  our  fitting  problem  with  Equation  (6)  the  amplitudes  go  in 
linearly,  but  not  so  the  geometric  parameters,  HX  and  SC.  So, 
we  have  to  face  a  non-linear  problem.  There  exist  different  me¬ 
thods  of  optimization  with  non-linear  parameters,  all  of  which 
are  necessarily  iterative.  The  criterion  expression  -  for  us: 
the  reduced  error  sum  -  as  function  of  a  number,  L,  of  diffe¬ 
rent  parameters  defines  an  L-dimensional  surface  in  an  (L+1 )- 
dimensional  space.  An  optimum  means  that  there  exists  a  minimum 
of  the  criterion  at  some  combination  of  parameters.  When  the 
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entries  are  eapirical  data,  there  is  not  necessarily  a  unique 
miniaua.  If  aore  than  one  appear,  we  encounter  a  serious  dif¬ 
ficulty. 

b.  Interdependence  of  parameters.  Another  problea 
which  aight  becoae  serious  occurs  when  the  paraaeters  are  not 
fully  independent  froa  each  other.  In  particular,  with  empiri¬ 
cal  data  such  connections  aight  only  be  valid  in  certain  ranges 
of  the  relevant  paraaeters.  In  that  case,  instead  of  a  clear 
and  unique  'crater-like'  minimum,  the  criterion  surface  shows 

a  aore  or  less  elongated  'minimum  trench'.  Note,  that  this  can 
in  practice  occur  without  strict  mathematical  interdependence; 
it  is  sufficient  that  there  exists  no  notable  change  of  the 
criterion  value  along  the  'trench'.  An  example  can  be  seen  on 
Figure  2  where  we  consider  only  one  member  in  the  sum  and  so 
have  L«2,  i.e.  only  two  variable  paraaeters,  HX  and  SC,  are 
adaitted.  Lines  of  constant  reduced  error  sum  S  in  an  HX  vs. 

SC  diagram  were  drawn.  We  call  such  diagram  'error  map'.  In 
such  case,  no  unique  optimization  result  can  be  reached.  The 
final  result  of  an  optiaization  procedure  can  then  be  any  point 
at  the  bottom  of  the  rift.  In  that  case,  the  final  result  de¬ 
pends  on  the  assumed  entry  values  of  the  iteration. 

c.  Admissible  number  of  unknowns.  The  iterative  pro¬ 
cedure  starts  usually  with  first  determining  the  direction  of 
steepest  descent  on  the  criterion- (error-)  surface  (arrows  in 
Figure  2).  L  direction  cosines  must  be  determined  to  this  end. 
If  the  assuaed  function  allows  it,  these  can  be  obtained  by 
functional  differentiation.  The  next  step  consists  in  a  first 
correction,  i.e.  changing  the  paraaeters  in  a  way  that  this  di¬ 
rection  is  consistently  followed.  A  step  width  SW  must  be  assu¬ 
med,  it  is  usually  decreased  with  increasing  order  of  the  ite¬ 
ration.  Though  this  looks  quite  promising,  the  execution  be¬ 
comes  aore  and  aore  involved  when  the  number  L  of  unknown  para¬ 
aeters  increases.  In  practice,  there  exists  a  limit  above  which 
the  iterative  process  might  be  unable  to  find  the  best  possible 
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set  of  parameter  values.  It  is  rarely  hopeful  to  let  L  go  up 
to  1o. 

In  our  case,  we  have  9  unknowns  if  the  number,  NFU,  of  LAX- 
functions  is  three,  12  if  it  is  4.  Further,  the  parameters  are 
not  independent.  In  particular,  the  amplitude  A  and  the  transi¬ 
tion  height  HX  are  de  facto  interrelated.  With  a  larger  A  and 
a  correspondingly  decreased  HX,  one  might  come  to  a  rather  si¬ 
milar  profile  as  with  the  original  values,  see  Figure  3.  Note, 
that  the  comparison  is  only  made  at  a  limited  number  (less  than 
M)  of  given  points  and  mainly  over  not  much  more  than  1  decade 
from  peak  value. 

d.  Failure  of  differentiating  procedure.  Therfore,  we 
could  not  apply  one  of  the  usual  optimisation  methods  unchanged. 
At  first,  we  tried  to  find  the  direction  of  steepest  descent, 
as  usual,  by  differentiation,  in  our  case  after  the  geometric 
parameters,  HX  and  SC.  When  doing  so,  the  functions  EPS0  and 
Eps^  are  needed.  Since  the  amplitude  A  -  noting  that  it  is  not 
independent  of  HX  -  was  held  constant,  the  so  found  HX  \SC  di¬ 
rection  was  not  really  that  of  steepest  descent  in  an  absolute 
sense.  In  fact,  by  adjusting  A  to  its  optimum  value  a  smaller 
error  sum  would  be  achieved.  On  the  other  hand,  due  to  the  in¬ 
terdependence,  a  method  admitting  all  three  variables  would  not 
lead  to  a  clear  answer,  because  the  iterative  process  might  de¬ 
velop  towards  extreme  values  of  one  or  the  other  parameters. 

With  this  result  in  mind,  we  finally  decided  to  build  an  opti¬ 
misation  procedure  of  our  own,  avoiding  functional  differen¬ 
tiation  and  taking  account  of  the  quasi-interdependence  stated 
above . 
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3.3  Our  Own  Optimisation  Procedure. 

a.  Elimination  of  the  amplitudes.  In  view  of  the  a- 
bove,  we  rearranged  the  iterative  process  applying  it  only  to 
the  geometric  parameters,  taking,  however,  care  that  the  opti¬ 
mum  amplitudes  for  all  three  or  four  LAY- functions.  Equation 
(6),  were  re-determined  at  each  step.  So,  independent  from  the 
applied  set  of  geometric  parameters,  the  amplitudes  were  step 
by  step  optimised  so  that  they  were  no  more  considered  as  vari¬ 
ables  in  the  iteration  itself.  Since  they  enter  linearly  in  F, 
these  coefficients  can  straightforward  be  found  by  resolving  a 
system  of  linear  equations.  With  three  or,  at  most,  four  un¬ 
knowns  this  needs  not  very  much  computing  time. 

b.  Parabolic  extrapolation.  Many  optimisation  pro¬ 
grams  apply  a  linear  extrapolation  along  the  line  of  steepest 
descent.  We  had  not  much  success  with  such  a  method,  because  it 
has  a  dangerous  tendency:  just  when  the  minimum  is  approached, 
the  extrapolation  towards  sero  error  pushes  the  iteration  to 
often  larger  aberrations  on  the  other  side  of  the  rift.  We, 
therfore,  finally  preferred  a  parabolic  extrapolation  towards 
sero  error . 

In  order  to  determine  this  parabola  (with  distance  on  the  line 
of  steepest  descent  as  abscissa),  we  need  three  points,  one 
more  than  for  a  linear  extrapolation.  We  divide  the  given  step- 
width  SW  by  2  and  make  two  consecutive  steps  along  the  given 
direction.  Together  with  the  start  position,  this  triple  deter¬ 
mines  a  parabola  which  either  is  concave  and  has  a  minimum 
somewhere  (inbetween  the  three  points  or  outside)  or  it  is  con¬ 
vex.  In  the  latter  case,  the  program  identifies  that  of  the 
three  positions  with  error  sum  lowest  and  then  takes  one 
more  step  (of  SW/2)  in  the  direction  beyond  this  point.  The  so 
found  fourth  position  together  with  the  two  others  next  to  it, 
is  then  used  as  another  triple  and  another  parabolic  extrapola¬ 
tion  is  undertaken.  The  same  is  done  when,  in  a  concave  situa¬ 
tion,  the  computed  minimum  is  too  far  distant  (situation  'FLAT': 
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for  us:  beyond  six  forward  or  two  backward  steps).  The  forward 
or  backward  shift  of  the  triple  is  arranged  by  the  subroutines 
SHIFTA  and  SHIFTR  ,  respectively.  The  number  of  these  shift  o- 
perations  (during  one  iteration  act)  must,  of  course,  be  limi¬ 
ted  (we  used  a  limit  of  5). 

Once  a  suitable  parabolic  minimum  is  found,  this  position  is 
taken  as  starting  point  for  the  next  iteration  act  (with  halved 
step  width). 

c.  Memorizing  best  conditions.  Since  it  is  not  cer¬ 
tain  that  this  position  has  really  a  better  error  sum  than  racked 
formerly,  during  the  whole  iteration  process  we  always  had  to 
remember  the  best  condition  found.  In  the  subroutine  RERRSS, 
which  computes  the  reduced  error  sum,  we  have  a  final  disjunc¬ 
tion  asking  whether  the  new  error  sum  is  smaller  than  the  smal¬ 
lest  one  reached  before  (which  is  in  the  memory).  If  it  is, then 
instead  of  the  memorized  set,  the  new  reduced  error  sum  toge¬ 
ther  with  the  corresponding  set  of  all  geometric  parameter  and 
amplitudes  is  commemorated.  In  case  no  improvement  was  reached 
the  next  iteration  act  starts  from  this  "provisionally  best" 
position  (with  renewed  determination  of  the  steepest  descent 
direction,  and  halved  step  width). 

d.  Finding  the  steepest  descent.  When  applying  diffe¬ 
rentiation  after  HX  and  SC,  we  took  the  amplitude  as  fixed  which 
is  not  justified.  This  could  be  corrected  for,  if  the  relation 
between  A^  and  HX^  a  priori  was  known.  But  it  is  not  and  its 
coefficients  depend  on  the  individual  data  set.  Therefore,  the 
direction  we  found  was  not  really  that  of  steepest  descent  when 
A  is  allowed  to  vary  too.  This  direction  error  was  found  to  be 
quite  appreciable  so  that  the  optimization  advanced  quite  slow¬ 
ly. 

In  order  to  get  this  direction  correctly  without  admitting  A  as 
an  independent  variable,  we  decided  to  use  a  non-differentiating 
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approach,  replacing  functional  differentiation  by  a  numerical 
one,  i.e.  by  determining  quotients  of  differences.  This  must, 
of  course,  be  done  with  adequate  small  steps.  At  start,  we  ta¬ 
ke  o. 5  and  1  km  (for  SC  and  HX,  respectively);  in  later  itera¬ 
tion  steps,  when  the  iteration  step  SW  decreases,  we  take  care 
that  the  'differentiation’  step  remains  small  against  it.  De¬ 
signating  the  parameter  values  by  Pj  (3  ;  1  ...  L,  according 
to  the  number  of  variables  chosen),  we  have  a  set  Pj  at  the 
start  position.  We  then  make  a  displacement  by  APj  in  just 
one  direction  j  so  that  we  get  a  set  of  L  auxiliary  points, 
each  on  one  of  the  ^-coordinate  axes. 

To  each  of  these  L  positions  the  error  sum  is  computed  (with 
optimised  amplitudes  always);  from  its  difference  against  that 
of  the  start  position  the  direction  cosine  (quotient  of  diffe¬ 
rences)  can  easily  be  computed.  The  £Pj  define  a  multidimen¬ 
sional  plane  and  its  direction  of  steepest  ascent  is  determi¬ 
ned.  This  does  subroutine  ASCENT.  The  direction  of  steepest 
descent  is  then  obtained  by  inverting  the  signs. 


e.  Description  of  the  main  program.  The  most  impor¬ 
tant  particularities  of  our  optimization  routine  APTR  have  been 
described  in  the  foregoing.  A  flux  diagram,  showing  the  most  im¬ 
portant  actions  in  this  program,  can  be  found  in  Figure  4a.  The 
important  subroutines  in  APTR  are  shown  in  the  subsequent  Figu¬ 
res  4.  The  full  text  of  subroutine  APTR  is  reproduced  as  Figure 
5  (in  computer  language  SIMULA, of  the  ALGOL  family). 


By  virtue  of  the  subroutine  CODELL,  different  combinations  of 
variables  can  be  selected.  Instead  of  applying  the  optimization 
routine  at  once  to  all  geometric  parameters,  it  was  in  most  con¬ 
ditions  more  advantageous  to  apply  it  several  times  with  each 
time  another  suitably  chosen  set  of  variables.  The  last  time 
we  then  admitted  all  geometric  parameters  as  variables,  but  now 
starting  from  a  rather  good  position  and  applying  a  small  step 
width.  In  case  of  NFU«3,  this  means  6,  for  NFU»4,  however,  8 
variables.  CODELL  allows  to  go  up  to  NFU-5  and  offers  2o  diffe¬ 
rent  combinations. 
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4.  RESULTS  MIDDLE  IONOSPHERE 

Our  data  source  for  this  height  range  are  W.  Becker's  profiles 
(see  chapter  2.1).  In  the  two  tapes  we  hare  received  from  this 
author,  for  each  profile  he  has  evaluated  he  gives  several  pre¬ 
sentations,  the  nos t  useful  for  out  purpose  being  one  in  height 
and  noroalized  plasma  frequency,  the  latter  in  units  of  1#  of 
the  critical  frequency  foP2,  thus  with  loo  points  per  profile. 
For  the  purpose  of  our  work  and  in  view  of  the  increase  of  com- 
puting  time  with  increasing  nunber  of  data  points,  we  retain 
only  every  fourth  data  point  which  makes  25  per  profile.  This 
should  be  enough  for  defining  a  profile  with  good  enough  accu¬ 
racy.  It  must,  however,  be  noted  that  even  so  the  lowest  samp- 
ling  frequency  (4%  of  foF2)  is  quite  a  low  one  on  which  actual 
ionosonde  measurements  could  actually  not  be  made:  either  -  at 
night  -  by  the  lower  frequency  limit  of  the  instrument,  or  -  by 
day  -  by  absorption  in  the  ionosphere.  In  fact,  Becker  has  ap¬ 
plied  extrapolation  towards  the  lowest  frequencies.  In  particu¬ 
lar,  at  night,  this  is  misleading  because  the  'nite  E'  layer  is 
known  to  exist,  but  does  not  show  up  on  the  ionogram  and  so  not 
in  Becker's  prof ilex.  Therefore,  we  have  given  zero  weight  to 
all  data  points  at  plasma  frequencies  below  o.5  MHz.  On  the 
other  hand,  a  decreased  weight  has  been  attributed  to  data 
points  at  altitudes  below  1o5  km  by  putting  W«Eps0(z; 1o5,4). 

As  explained  in  our  earlier  report,  there  are,  amongst  Becker's 
profiles,  quite  clearly  cases  where  something  went  completely 
wrong  with  the  height  scale  or  with  extrapolation  so,  that  geo¬ 
physically  improbable  positions  of  the  E-layer  appear  from  time 
to  time.  Or,  this  layer  is  known  to  be  the  most  regular  of  all 
and  to  appear  in  the  same  height  range  always.  These  errors  are 
probably  due  to  ionogram  scaling  or  occurred  during  digitizing. 

4.1  Night  Profiles.  At  night,  lonosondes  only  see  an  F-regi- 
on  trace  with  a  minimum  true  height  well  above  loo  km.  It  is, 
however,  known  from  other  sources  that  an  E-reglon  is  always 


existing,  but  with  rather  low  peak  electron  density  so,  that 
it  cannot  be  readily  detected  with  a  normal  ionosonde.  Piggott 
and  Rawer1  give  values  of  foE  which  lie  between  l.o  and  o.55 
MHz  /Ne  between  12.4  and  3.75  •  1o9  n This  should  held 
from  1  h  after  sunset  to  Wmi before  sunrise.  Becker's  tables, 
however,  give  curves  ending  with  1#  of  foF2  at  heights  well 
above  loo  km,  the  lowest  points  of  which  have  certainly  been 
obtained  by  extrapolation  and  cannot  be  correct.  In  our  compu¬ 
tations  we  always  disregard  those  points  of  Becker's  at  which 
the  plasma  frequency  is  below  o.5  MHz. 

a.  Odd  conditions.  At  tines,  the  bottom  part  of  a 
profile  nay  show  a  steep  gradient,  at  other  tines,  however,  an 
almost  linear  decrease  of  plasna  frequency  with  height.  We  have 
sone  doubts  about  these  latter  cases,  because  it  is  known  that 
a  deep  valley  exists  at  night  so,  that  a  sharp  gradient  should 
appear  sonewhere  in  the  2'2-profile.  Our  doubts  are  increased  by 
the  observations  that  such  cases  often  appear  for  a  short  peri- 
or  anongst  'regular'  ones.  Looking  through  a  dense  series  of 
sunner  night  measurements  (about  every  half  hour),  we  found  ma¬ 
ny  cases  of  "switching"  between  the  two  types  in  a  3o  min  in¬ 
terval.  Apparently,  some  non-stationary  condition  of  the  F- 
region  must  have  occurred.  May-be,  as  by  day,  gravity  waves  pro¬ 
duce  short-lived  deformations  of  the  ionogram  trace  which,  at 
data  reduction,  are  smoothed  and  appear  in  the  final  profile  as 
an  extension  towards  lower  altitudes. 

This  odd  kind  of  night  profiles  in  Becker's  set  has  a  large 
range  of  almost  linear  variation  of  plasma  frequency  with 
height.  Figure  6  shows  a  normal  and  an  odd  profile  for  compari¬ 
son. 

1  W.R.  Piggott  and  K.  Rawer,  URSI  Handbook  of  Ionogram  Inter¬ 
pretation  and  Reduction,  Elsevier,  Amsterdam  1961,  p.  127. 


When  fitting  with  Just  one  LAY-function,  the  normal  profiles 
use  to  be  very  well  reproduced,  the  odd  ones  not  so  well.  This 
appears  clearly  from  the  corresponding  error  maps  which  we  show 
in  Figure  7.  The  normal  profile  of  Figure  6a  shows  a  clear-cut 
'trench*  in  the  error  map  (Figure  7a).  Since  HX  and  A  are  in¬ 
terrelated,  the  optimum  condition  can  be  realized  with  diffe¬ 
rent  combinations  of  both  as  shown  by  the  'trench'.  The  odd 
profiles  (Figure  6b),  however,  show  instead  of  a  'trench'  a 
broader  'valley'.  The  absolute  scale  values,  SC,  are  much  lar¬ 
ger  in  odd  cases,  more  than  twice  those  found  with  parabolic 
shapes.  This  stems  from  the  fact  that  our  "best  fit"  is  one 
which  covers  the  over-all  width  of  the  layer  rather  than  to 
reproduce  best  its  high-density  (peak)  range. 

As  expressed  above,  we  have  doubts  whether  these  conditions 
which  occur  only  for  short  time  sections  should  be  taken  as 
representative.  Rather,  we  shall  Dase  our  summary  on  "quasi- 
parabolic"  conditions  only.  This  is  the  more  Justified  since 
we  should  anyway  disregard  the  doubtful  information,  given  at 
xower  frequencies  in  Becker's  profile.  In  the  synthesis  of 
night-time  middle  ionosphere  profiles,  we  might  then  introduce 
a  deep  valley  and  "nite-E". 

Just  in  order  to  demonstrate  how  the  normal  and  the  odd  types 
behave,  in  Figures  7  we  have  shown  error  maps  (arbitrary  units) 
for  these  two  conditions.  Both  are  only  3o  min  apart.  The  cor¬ 
responding  profile  shapes  can  be  seen  from  Figures  6.  The  pro¬ 
file  at  2oh3o  is  odd,  that  at  21  h  normal.  The  quick  change  of 
types  might  be  seen  from  the  following  table: 
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Table  1:  Variation  of  scale,  SC,  for  two  night  periods 
(25/27  June,  1954).  NFU  -  Is  one  LAY-function 
only.  (Clearly  odd  profiles  underlined). 


(a)  Night  25/26  June: 


Hour 

19.5o 

2o 

2o.3o 

21 

23 

23.3o 

00.3o 

01 

SC/km 

23 

14 

2 4 

8 

4 

6 

23 

2o 

(b)  Night  26/27  June: 


Hour 

2o.27 

2o.57 

22 

23.57 

00.27 

SC/km 

4o 

24 

1o 

11 

21.5 

It  is  a  merit  of  Becker's  dense  series  of  profiles  that  such 
phenomena  and  their  changes  can  be  seen  at  all  which  does, 
of  course,  not  hold  for  averaged  profiles.  We  feel  that,  in 
view  of  our  final  aim,  we  have  to  disregard  these  cases  of 
'linear  decrease'  because,  what  counts  for  our  purpose,  are 
stable  conditions. 

On  the  other  hand,  comparing  with  density-averaged  profiles, 
the  individual  profiles  in  Becker's  set  have  somewhat  smaller 
thickness.  We  conclude  that  averaging  over  all  conditions  leads 
to  thicker  layers  than  does  a  selection  of  stable  profiles. 
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b.  Normal  conditions  at  night.  Apart  from  short-lived 
deformations,  normal  night  profiles  use  to  be  of  almost  ’ para¬ 
bolic'  shape.  This  is  valid  at  all  latitudes.  Therefore,  for 
reference  purpose  at  least  the  night  profiles  can  be  reproduced 
by  just  one  LAY- function.  When  trying  with  two  such  functions, 
one  amplitude  happened  to  come  out  quite  small  so  that  the  ef¬ 
fect  of  the  second  member  in  Equation  (6)  was  negligible.  This 
might  be  different  with  odd  shapes  which  sometimes  do  occur 
(see  a.  above).  We,  therefore,  feel  definitively  that  the  night¬ 
time  F-region  profile  should  be  reproduced  with  just  one  such 
function. 

One  LAY-function  has  only  two  geometric  parameters,  such  that 
our  mapping  program  gives  the  best  answer  directly  without  fur¬ 
ther  assumptions.  For  the  majority  of  profiles  which  are  of  al¬ 
most  parabolic  shape,  the  HX  versus  SC  map  shows  a  straight 
"trench"  of  "valley"  which  lies  vertically  in  the  map,  such 
that  the  scale  where  the  minimum  appears  is  well-defined.  In  the 
trench  one  obtains  about  the  same  error  value  with  different 
transition  heights  HX.  This  structure  is  the  simplest  we  found. 
At  analysis,  this  means  that  while  the  scale  must  be  optimised 
the  transition  height  is  left  free  in  a  large  range  (of  loo  km 
or  more). 

This  feature  appears,  of  course,  only  with  our  schedule  of  re¬ 
determination  of  the  amplitude  A  at  each  step.  With  a  fixed  am¬ 
plitude  one  would  obtain  just  one  minimum,  not  a  trench.  Some 
examples  of  error  maps  can  be  found  in  Figures  8  (which  we  shall 
discuss  below). 

Along  the  trench  there  exists,  of  course,  a  quite  appreciable 
variation  of  the  amplitude  A  which  decreases  towards  larger  va¬ 
lues  of  HX  (downwards  in  the  Figures).  The  'trench*  is,  however, 
not  unlimited  in  the  HX-direction :  at  its  lower  end  (large  HX) 
there  appears  a  quick  increase  of  the  error  sum.  At  the  other 
end,  the  amplitude  becomes  larger  and  larger  and  at  a  certain 
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limit  there  appears  an  even  steeper  increase  of  the  error  sub. 
So,  for  a  given  profile  we  have  to  accept  a  limited  HX  range 
with  a  rather  well  defined  scale  SC  as  optimum  condition.  In 
practice,  it  is  preferable  to  have  not  too  large  amplitudes 
and,  therefore,  we  shall  propose  for  applications  HX  values 
near  the  lower  end  of  the  trench. 

All  night-time  error  maps  (of  which  we  have  drawn  a  large  num¬ 
ber)  show  similar  structure,  though  the  width  of  the  trench  va¬ 
ries.  It  is  important  to  note  that  the  low  amplitude  limit  (at 
the  bottomside  in  our  Figures)  was  always  found  well  below  the 
peak  altitude  HM,  sometimes  even  below  HM/2. 

As  for  our  optimisation  program,  it  reaches  necessarily  a  cor¬ 
rect  and  unique  scale  value  SC,  while  the  final  value  of  HX, 
due  to  the  particular  configuration,  usually  remains  near  the 
start  value  of  this  parameter.  This  ambiguity  is  typical  at 
profile  analysis. 

At  profile  synthesis,  however,  (which  is  later  to  be  executed 
when  applying  our  results  in  IRI)  we  have  quite  a  bit  of  free¬ 
dom  in  choosing  HX  but  not  so  for  SC.  Since  this  is  our  final 
intention,  we  have  studied  in  more  detail  where  the  low  ampli¬ 
tude  end  of  the  trench  occurs  and  whether  the  relevant  HX-value 
can  be  specified.  To  this  end,  that  value  (designated  as  HX1 ) 
was  compared  in  turn  with  the  peak  altitude,  HM,  the  half-den- 
sity  altitude,  HHA,  and  the  quarter  density  altitude,  HQA.  In 
Figures  9a,  b,  c  the  relevant  ratios  have  been  plotted  against 
the  peak  altitude,  HM.  Data  were  from  the  solar  minimum  year 
1954  on  the  one  hand,  and  from  the  very  high  solar  maximum  in 
1957  and  1958  on  the  other.  The  Figures  show  the  ratios  decrea¬ 
sing  with  increasing  HM.  A  straight  line  connection  of  the  me¬ 
dian  point  of  the  two  periods  is  given  by  the  following  Equa¬ 
tions  : 


(♦):  HX^/HM  -  .795o  -  ,59o  HM/looo 


(6a) 


3o 


(o): 

HX^/HHX  -  .981o  - 

.69o  HM/looo  km 

(8b) 

(x): 

HX-j/HQU  -  .96o5  - 

.525  HM/looo  km 

(8c) 

The  dispersion  of  the  individual  points  is  rather  large.  Never¬ 
theless,  the  decreasing  tendency  is  certain.  We  have  also  plot¬ 
ted  the  corresponding  height  differences  (in  Figures  loa,  b,  c). 
This  is ,  of  course,  increasing  with  increasing  HM  and  the  con¬ 
nection  of  the  median  points  is  given  by: 


(  +  ): 

HM  -  HX1  «  .6130 

HM  -  7o.3o 

km 

(9a) 

(o): 

HHA  -  HX1 «  .3680 

o 

• 

00 

■V 

1 

km 

(9b) 

(x): 

HQU  -  HX.j»  .2845 

HM  -  38.65 

km 

(9c) 

We  leave  it  open  which  of  these  relations  might  be  most  helpful 
at  profile  synthesis.  The  above  relations  (..a)  with  HM  have 
the  advantage  that  they  do  not  need  another  characteristic. 

HQU  might  not  be  available  experimentally  when  the  peak  density 
is  small.  So,  the  choice  would  probably  be  HM  (or  HHA).  However, 
since  these  are  median  lines  half  of  our  data  ly  on  the  unwanted 
side.  For  this  reason,  at  applications  we  should  not  use  the  me¬ 
dian  line  but  one  on  the  safer  side,  e.g.  the  appropriate  quar- 
tile  line  (thin  lines  in  the  Figures).  These  obey  the  following 
Equations: 


Hfy  /HM 
i.e.  HXq1 

HM-HXq2 

i.e.  HXq2 


0.7184  -  o.57o8  HM/looo  km, 
0.7184  *HM  -  o.57o8  HM2/1ooo  km 

0.6  HM  -  So  km 

0.;<  HM  +  30  km 


} 

1 


(loa) 


(lob) 


Since  it  has  been  proposed  to  apply  a  simple  proportionality  re- 


-  31 


lation  for  deducing1  HHA  from  HM,  we  investigated  this  ratio 
with  our  data.  Figure  11  shows  the  result.  For  night  conditions, 
in  solar  minimum  and  maximum  as  well f the  ratio  is  quite  near  to 
0.8  as  indicated  by  this  author.  Drawing,  as  above,  the  connec¬ 
ting  line  between  the  median  points  of  both  periods,  we  get 

HHA/HM  -0. 9o75  -  0.255  HM/looo  km. 

Since  the  dispersion  is  rather  small,  an  oblique  median  line 
through  all  points  should,  however,  be  preferred;  it  reads: 

HHA/HM  -0.91  -0.25  HM/looo  km  (11) 

This  is  the  bold  line  in  Figure  11.  The  quartile  ranges,  as  in¬ 
dicated  by  thin  lines  in  the  same  Figure,  show  that  the  depen¬ 
dence  on  HX  is  significant.  During  solar  minimum  when  the  peak 
altitude  used  to  be  considerably  smaller,  the  median  value  of 
the  ratio  was  o.83  while  only  0.785  during  solar  maximum. 

c.  Short  term  variations.  Becker’s  data  set  covers 
some  dates  for  which  profiles  were  computed  at  hourly  or  even 
half-hourly  intervals.  These  are  well  suited  for  answering  the 
question  after  the  short  term  variability  of  the  parameters  HX 
and,  in  particular,  SC.  In  Figure  12  we  show  in  some  detail  what 
happened  during  the  late  night  of  24/25  June,  1954.  Apart  from 
the  fp  (plasma  frequency)  vs.  height  profile,  two  cuts  of  the 
error  map  (around  the  minimum)  are  shown  with  HX  fixed  at  HM/2 
and  HM/3,  respectively.  It  appears  from  this  Figure  that  the 
quick  changes  between  "normal”  (a,c,f,g,h)  and  "odd"  (b,d,e) 
profile  types  are  accompanied  in  the  map  by  a  switching  between 
a  narrow  "trench"  centered  at  a  small  SC  (between  5  and  25  km) 
and  a  broader  one  with  much  larger  SC. 


T.L.  Gulyaeva,  Implementation  of  a  new  characteristic  para¬ 
meter  into  the  IRI  sub-peak  electron  density  profile, 

Adv.  Space  Res.  2,  No.  1o,  191-194  (1983). 


One  night  try  to  take  care  of  the  different  behaviour  of  what 
we  called  above  "normal"  and  "odd"  layers  by  comparing  with  a 
layer  thickness  parameter,  e.g.  the  thickness  SH,  given  by  HHA 
(the  point  of  half  peak  density),  or  SQ,  given  by  HQU  (at  half 
plasma  frequency,  i.e.  quarter  density).  These  points  are  shown 
on  the  profiles  of  Figure  12  (as  cross  and  bar,  respectively). 

It  appears  that  SH  »  HK  -  HHA  is  not  much  different  while  the 
quarter-thickness  SQ  ■  HM  -  HQU  is  considerably  larger  when¬ 
ever  the  layer  shape  is  "odd".  The  drawings  also  show  the  pre¬ 
sumed  "typical"  value  of  0.8  HM  for  Nm/2.  It  appears  that  -  at 
night  -  this  altitude  is  almost  always  below  HHA  such  that, 
when  applying  the  coefficient  0.8,  one  makes  the  layer  thicker 
than  observed.  /Compare  also  Figure  11  above/. 

We  have  studied  in  more  detail  another  night  of  the  same  month 
/ 26/27  June  1954/  producing  error  maps  and  cuts  at  HX1  «  HM/2 
for  15  profiles  which  Becker  has  deduced  for  the  period  l8h3o 
through  o5h3o  LT.  It  appeared  that,  when  arranging  after  the 
optimum  SC  value,  we  find  three  groups  of  profiles:  a  group  (a) 
with  SC  between  5  and  12  km,  a  group  (b)  around  2o  km  and  a 
group  (c)  above  25  km.  At  least  in  that  particular  night  the 
groups  are  well  distinguished.  We  show  the  error  map  cuts  in 
Figures  13a,  b,  c  according  to  these  groups.  (The  order  inside 
each  group  was  chosen  for  convenience  of  the  aspect  of  the  dra¬ 
wing).  Of  course,  the  trenches  are  narrower  in  group  (a)  where 
the  optimum  SC  is  small.  The  width  can  become  quite  large  in 
group  (c).  This  means  that  we  have  a  large  range  of  acceptable 
SC  values,  while  in  group  (a)  the  "optimum  analysis"  leads  to  a 
well  defined  value.  The  narrowest  trench  was  found  for  o2  h 
with  a  quite  low  minimum  error  which,  however,  increases  rapid¬ 
ly  at  both  sides. 

The  main  question  on  our  context  is  whether  on  behalf  of  the  op¬ 
timum,  SC(opt),  there  exists  some  regularity  or  a  systematic  va¬ 
riation  during  the  night.  In  order  to  demonstrate  that  this  is 
not  so,  we  show  in  Figure  14  some  characteristic  features  of  the 
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error  map  cuts  in  chronological  order.  On  the  diagram  on  top, 
showing  SC(opt),  one  easily  identifies  the  three  groups.  Their 
appearance  during  that  night  is  apparently  at  random,  no  syste¬ 
matic  variation  being  visible.  Also,  there  exists  no  relation 
whatsoever  with  the  HX-j  value  (HM/2,  shown  on  the  bottom  dia¬ 
gram).  The  same  holds  for  the  minimum  reduced  error  sum.  Err, 
which  was  plotted  in  second  position  (in  units  of  1o-*).  As 
mentionned  above,  there  must  exist  a  relation  with  the  width 
of  the  trench.  In  Figure  14  we  use  two  different  definitions  of 
width:  ^  is  the  width  for  twice  the  minimum  Err-value,  while 
£  2  is  width  at  the  fixed  value  Err«5o  •  1o“  .  For  future 
applications  at  profile  synthesizing,  <T 2  night  be  of  greater 
interest  because  it  means  considering  the  same  maximum  accep¬ 
table  error  sum  for  all  profiles  equally. 

For  synthesizing  application  one  might  yet  be  interested  to  have 
the  median  values.  In  the  particular  case  of  Figure  14  these  are 
found  as: 

HX1  «  15o  km;  SC(opt)  ■  21  km;  Width  *  9  km. 

Summarizing  our  findings  /in  view  of  future  synthesizing/  we 
may  state  that,  at  night,  types  and  characteristics  of  the  F2- 
profile  are  quickly  variable  in  an  irregular  manner.  Since,  ap¬ 
parently,  conditions  can  change  seriously  in  an  interval  of 
3o  min,  it  is  not  helpful  to  ask  for  a  regular  variation  with 
the  hour.  Therefore,  median  parameters  should  be  adopted  at 
night,  which  might,  of  course,  yet  depend  on  season  and  solar 
activity. 

In  this  context,  it  might  be  indicated  to  apply  a  HX1  value 
which  follows  the  peak  altitude  HM.  In  the  Figures  above  we 
simply  used  HM/2  (and  also  HM/3  in  Figures  12).  It  might,  how¬ 
ever,  be  preferable  to  use  the  relation 

HX1  -  HM  -  A 

where  ^  could  be  taken  from  Equation  (\  ob)  •  (During  the  solar 
minimum  year  1954,  this  value  was  rather  near  to  HM/2). 


d.  Long  term  variations.  We  have  quite  generally  found 
that,  at  night,  for  each  individual  profile  which  we  have  ana¬ 
lyzed  the  optimum  scale  SC (opt)  is  well  defined,  though  the 
statistical  dispersion  of  all  profiles  is  rather  large.  Going 
over  Becker's  profiles  for  each  night  of  his  data  set,  we  stu¬ 
died  a  few  around  midnight.  As  mentionned  above,  che  main  data 
are  from  summer  1954  and  from  spring  1958.  In  Figures  8  we  have 
shown  eight  typical  error  maps  selected  from  the  whole  set.  It 
appears  that,  though  the  position  of  the  optimum  SC-value  can 
be  largely  different,  the  structure  of  these  maps  is  always 
similar.  It  is  characterized  by  a  "trench"  extending  along  the 
ordinate  (HX)  direction,  but  the  extension  of  which  is  clearly 
limited  at  both  ends.  The  end  with  largest  HX  admissible  (bot¬ 
tom  side  in  the  Figures)  has  the  smaller  amplitudes  and  is  of 
particular  interest  for  future  synthesizing  of  profiles.  (This 
end  was  studied  in  Figures  9  and  1o,  Subsection  b). 

It  should  be  noted  that,  by  night,  our  weighting  covers  all  of 
Becker's  profile  points  down  to  a  plasma  frequency  of  o.5  MHz. 
It  is  evident  that  another  fit  of  the  profile,  limited  to  a 
small  height  ranqe  near  the  peak,  might  end-up  with  somewhat 
other  values.  Just  this  profile  range  is,  however,  in  Becker's 
analysis  determined  by  an  assumed  nose  function  (mostly  a  para¬ 
bola).  For  a  reference,  it  is  certainly  more  important  to  fit 
the  profile  as  whole,  what  we  have  done. 

In  order  to  recommend  an  SCI,  we  first  felt  that  SCI  could  be 
related  with  the  half-density  thickness  SH  (or  the  quarter-den¬ 
sity  thickness  SQ).  Figure  15  shows  correlograms  of  the  pairs 
SH  vs.  SC  and  SQ  vs.  SC.  There  is,  apparently,  some  correlation 
for  SC-values  below  about  4o  km.  For  higher  values,  however,  we 
cannot  Bee  any  correlation. 

It  might,  thus,  be  the  best  solution  to  recommend  the  median  SC 
of  the  monthly  sets.  Se,  we  find  (around  midnight): 

SCI  «  2o  in  summer  1954  (solar  minimum),  but 

SCI  »  45  in  spring  1958  (solar  maximum). 


It  mist  yet  be  found  out  whether  the  effect  of  the  solar  cycle 
or  that  of  the  season  is  more  important. 

One  might  find  it  more  appropriate  to  link  the  SCI -value  with 
the  peak  altitude  -  which  is  much  higher  during  solar  maximum. 
Using  the  HM  group  medians  from  Figures  9  and  1o  which  are  3 14 
and  481  km,  respectively,  and  assuming  a  linear  relation,  we 
find 

SCI  «  o.15  HM  -  27.1  km  (12) 

as  a  median  relation  to  be  used  by  night. 

As  for  HX1 ,  as  we  have  seen,  it  can  be  chosen  in  a  large  range 
but  must  not  come  near  to  HM.  In  Subsection  (b)  we  have  investi¬ 
gated  quite  a  few  error  maps  determining  the  greatest  admissible 
value  of  HX1 .  We  propose  to  apply  the  quartile  value  of /l  ,  as 
specified  in  Equation  (19b).  This  ends  up  with  a  recommended  HX1 
after 

HX1  »  0.<*  HM  +  30  km.  (13) 

This  relation  gives  HX1  «  155  a®  the  1954  (summer)  median  and 
223  km  for  spring  of  1958,  the  ratio  to  HM  being  o.$ro  and  o.^g, 
respectively. 

It  should  be  noted  that  the  above  two  standard  parameters  do  not 
yet  fully  determine  a  profile.  One  more  parameter  is  left,  name¬ 
ly  the  amplitude  A.  At  synthesis,  it  could  be  chosen  so  as  to 
satisfy  one  constraint,  e.g.  the  half  or  quarter  density  point. 
Since  the  peak  of  &  LAY-function  is  fixed,  the  profile  goes 
through  a  predetermined  HHA  or  HQU,  provided  A  is  chosen  by  the 
condition 

A  -  -o.3o1/LAY  (HHA;  HX,SC)  (14a) 

or  A  -  -o.6o2/LAY  (HQU;  HX,SC).  (14b) 

It  is,  of  course,  important  to  find  geophysical  relations  de¬ 
termining  HHA  and  HQU. 


A  future  night  reference  must,  of  course,  also  specify  a  val¬ 
ley  and  an  E-region.  These  do  not  appear  on  ionograms.  The  re¬ 
levant  information  must  be  taken  from  other  sources  (as  in  the 
present  IRI).  From  an  extended  study  of  night  ionograms  (mid¬ 
latitude)  Piggott^  established  the  following  table  for  the  ave¬ 
raged  night  values  of  foE,  depending  on  the  time  since  sunset 
or  before  sunrise.  One  starts  from  an  observed  foE  value, 
reached  near  sunset  (first  column  in  Table  2a),  and  goes  to¬ 
wards  a  corresponding  value,  reached  shortly  after  sunrise 
(last  column  in  Table  2b). 

Table  2.  Night  Values  of  E- Layer  Critical  Frequency/KHz 
(after  Piggott). 


a.  h  =  hours  after  sunset  (with  first  value  of  foE) 


h 

0 

o.25 

o.5 

o.75 

1  .0 

1.5 

2.o 

2.5 

3.o 

3.5 

4  .o 

o 

• 

A 

1  .6 

1.35 

1 .2 

1.1 

1  .0 

0.85 

o  .8 

o.75 

o.7 

o.65 

o.6 

o.55 

1.5 

1 .3 

1.15 

1  .o5 

1  .0 

o.85 

o.8 

o.75 

o.7 

o.65 

o.6 

o.55 

foE 

1.4 

1.25 

1.1 

1  .0 

o.95 

o.85 

o.75 

o.7 

o.65 

o.6 

o.6 

o.55 

1.3 

1.15 

1  .o5 

1  .0 

o.9 

0.8 

o.75 

o.7 

o.65 

o.6 

o.6 

o.55 

1.2 

1.1 

1  «o 

o.9 

o.85 

o.8 

o.75 

o.7 

o.65 

o.6 

o.55 

o.55 

1.1 

1  .0 

o.95 

o.9 

o.85 

o.75 

o.7 

o.65 

o.6 

o.6 

o.55 

o.55 

b.  h  =  hours  before  last  value 


h 

>2.o 

2.o 

1.5 

1  .0 

o.75 

o.5 

o.25 

C 

o.55 

o.65 

o.75 

o.9 

o.95 

1.1 

1 .3 

1.6 

o.55 

o.65 

o.75 

0.85 

o.95 

1  .o5 

1  .2 

1.5 

foE 

o.55 

o.65 

o.75 

o.85 

o.9 

1 .0 

1.15 

1 .4 

o.55 

o.65 

o.75 

0.8 

o.9 

o.95 

1.1 

1.3 

o.55 

o.65 

o.7 

0.8 

o.85 

o.95 

1  .o5 

1.2 

o.55 

o.65 

o.7 

o.75 

0.8 

o.9 

o.95 

1.1 

W.h.  Piggott  and  K.  Rawer,  URS1  Handbook  of  Ionogram  Inter¬ 
pretation  and  Reduction,  Elsevier,  Amsterdam  1961. 
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Determining  the  night-E  peak  by  the  values  given  in  Table  2  and 
using  standard  values  for  the  peak  height  hmE  (1o5  km)  and  for 
the  valley  bottom  (e.g.  14o  km),  we  yet  need  a  valley  depth, 
e.g.  log(NmV/NmE) .  It  is  known  that  the  night  time  E-F-valley 

is  quite  deep,  but  the  estimates  cover  a  large  range.  In  IRI  the 

2 

average  results  of  Maeda's  study  are  used  to  this  end,  though 
they  are  based  on  rather  old  rocket  measurements  the  calibration 
of  which  might  be  doubtful.  This  is  still  a  difficult  problem. 

More  recent  investigations  with  the  coherent  scatter  sounding 
technique  are  reported  to  have  shown  that,  at  night,  in  the 
broad  valley  range  rapid  and  important  short-term  changes  occur 
quite  regularly.  Transient  (mostly  descending)  irregular  layers 
might  reach  considerably  increased  electron  density  than  ave¬ 
rage.  For  a  reference,  one  should  disregard  such  features  and 
better  give  a  quasi-minimum  value. 

Once  the  NmV/NmE  ratio  is  fixed  and  standard  geometric  parame¬ 
ters  are  assumed  for  the  second  and  third  LAY- function,  the 
night-F  profile  might  be  completed  down  to  hmE.  Since  the  main 
LAY-function  represents  the  observable  F2-layer  down  to  low 
density,  the  two  other  functions  can  be  determined  with  fixed 
geometric  parameters,  the  two  amplitudes  being  determined  by 
three  or  four  conditions (constraints).  A  solution  can  be  found 
by  applying  the  method  of  least  squares.  This  is  outside  of  the 
present  study. 


2 

K.-I.  Maeda,  Study  on  electron  density  profile  in  the  lower 
ionosphere,  J.  Geomagn.  Geoelectr.  23,  133-159  (1971). 
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4.2.  Day  profiles.  The  fitting  situation  is  quite  dif¬ 
ferent  with  day  profiles.  At  night,  the  E-region  is  so  poorly 
ionized  that  the  P-region  could  independently  be  fitted.  This 
was  due  to  the  fact  that  there  was  no  overlap  with  the  range 
of  the  main  LAY-function.  By  day,  however,  for  an  optimum  fit 
the  ranges  of  the  different  LAY-functions  are  alnost  always 
overlapping.  Therefore,  the  P-region  profile  cannot  independent¬ 
ly  be  fitted;  the  fitting  procedure  has  to  take  account  of  all 
three  or  four  functions  and  their  parameters.  This  means  that 
we  have  six  or  eight  geometric  parameters  which  must  be  opti¬ 
mized.  The  methods  we  applied  to  this  end  were  described  in 
Chapter  3. 

There  occur  two  classes  of  day-time  profiles  which  must  be  dis¬ 
cussed:  when  no  Pi-layer  is  visible,  the  middle  ionosphere  pro¬ 
file  can  quite  well  be  fitted  with  three  LAY-functions.  The 
main  function  which  has  the  largest  scale  SCI  reproduces  the 
main  (P2-)  peak  but  covers  the  whole  range.  The  secondary  func¬ 
tions  are  needed  to  (i)  reproduce  a  valley  (or  a  turning  point) 
somewhere  around  14o  km,  and  (ii)  the  peak  of  the  E-region  at 
about  1o5  km.  These  have  smaller  scales  but,  since  their  HX  va¬ 
lues  are  rather  near  together,  they  have  overlap  with  each 
other  (and,  of  course,  with  the  main  function). 

The  appearance  of  the  PI -feature  changes  the  situation.  The 
corresponding  turning  point  in  the  F1/P2  profile  can  only  be 
reproduced  with  a  fourth  LAY-function. 

a.  Profiles  without  the  PI  feature.  These  are  represen¬ 
ted  by  one  main  layer  (P)  and  two  secondary  ones  (V  and  E), 
thus  with  I  ■  3  /in  Equation  (6)/.  Pigures  16a, b  show  typical 
examples  where  an  FI -layer  did  not  appear  at  all.  Figure  16c, 
however,  shows  FI,  but  this  feature  -  by  suitable  starting  con¬ 
ditions  -  was  neglected  at  fitting.  Apart  from  the  height  range 
where  PI  occurs,  the  profile  is  certainly  representative.  A 
larger  selection  of  profiles  fitted  with  I  »  3  is  reproduced  in 
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Figures  17.  It  cannot  be  denied  that  these  profiles  are  all  re¬ 
produced  quite  satisfy ingly.  So,  1=3  is  the  right  choice  un¬ 
der  the  assumed  conditions.  It  is  also  apparent  that  the  HX  and 
SC  values  of  the  second  and  third  LAY-function  are  quite  stable. 

The  median  geometric  parameters  of  these  secondary  layers  were 
found  from  a  statistical  evaluation  of  a  number  of  day-time 
profiles  obtained  in  spring  of  1958  (a  year  of  very  high  solar 
activity).  The  median  values  found  for  this  period  are: 

HX2  SC2  HX3  SC3 

116  2  1o5  5  km.  (i) 

As  for  the  parameters  of  the  main  layer  (1),  there  exists  some 
interpendence  as  we  have  mentionned  earlier.  Thus,  there  is  a 
certain  freedom  in  choosing  HX1 ;  in  the  fitting  procedure  the 
result  depends  largely  on  the  starting  value.  It  appeared  rea¬ 
sonable  to  admit  proportionality  with  the  peak  altitude  HM.  We 
used  to  start  fitting  with  a  ratio  HX  /HM  =  o.9,  and  this  was 
essentially  conserved  during  the  fitting  procedure.  So,  we  have 
with 

HX1  o.9  HM  and  o.85  HM  (ii) 

found  the  median  thickness  scale  to  be 

5CT  «  63  km  and  7o  km,  respectively.  (iii) 

Once  the  secondary  LAY-parameters  have  been  optimized  /as  in 
(i)/,  one  may  proceed  to  draw  an  error  map  (with  H  vs.  SCI  as 
coordinates).  The  details  of  the  map  are,  of  course,  somewhat 
depending  on  the  parameters  taken  for  the  secondary  functions. 
There  are,  however,  two  typical  features  occurring  rather  regu¬ 
larly:  one  is  that  -  quite  different  from  night  conditions  - 
there  ixists  now  a  lower  limit  for  HX1  (on  top  in  our  figures); 
secondly,  there  is  again  a  trench  which,  however,  extends  up  to 
HX1  values  far  beyond  HM.  The  lowest  Err-values  are  now  found 
when  HX  is  rather  near  to  HM.  Examples  are  shown  in  Figure  18. 
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As  done  above  for  night  (Figure  15),  we  established  correlo- 
grans  between  SCI  and  the  half-density- thickness  SH  and  quarter 
density  thickness  SQ.  For  the  considered  period  median  values 
and  dispersion  (decile)  ranges  were  found  as  shown  in  Table  3* 
Looking  for  correlations,  none  could  be  seen  with  SQ.  With  SH, 
though  a  median  regression  line  was  empirically  found  as 
SCI  *  5  SH  -  525  km,  it  could  not  be  taken  as  a  satisfactory 
description  since  it  gives  negative  values  of  SCI  for  SH  below 
1o5  km.  Although  a  line  from  the  origin  to  the  center  of  the 
distribution  around  this  line  is  not  symmetric  but  oblique. 

With  other  words,  the  ratio  SC1/SH  is  not  constant  but  decrea¬ 
ses  with  increasing  value  of  SH.  We  felt  that  a  more  involved 
relation  had  to  be  established  taking  account  of  this  feature. 
By  log  vs.  log  plotting  we  found: 

SCI  «  (SH/42  km)4 

This  formula  also  takes  account  of  the  fact  that  the  ranges  over 
which  both  parameters  vary  are  rather  different,  as  can  be  seen 
from  Table  3. 

Table  3.  Thickness  parameters  (F-region,  day,  1*3 ) 

SH/km  SQ/km  SCI /km 

median  123  187  63 

quartile  range  114-127  159-2o4  49-88 

decile  range  9o-15o  132-22o  3o-13o 


b.  Profiles  with  FI.  We  have  shown  that,  generally,  re 
production  of  a  day-profile  with  three  LAY-functions  (I»3)  is 
largely  satisfying.  We  have  shown  by  Figure  16c  that  this  is 
true  even  with  FI  present,  provided  one  accepts  some  deviation 
from  the  original  profile  in  a  rather  limited  height  range. 


However,  if  it  is  intended  to  reproduce  the  FI  feature  as  such, 
a  typical  summerday  profile  cannot  accurately  be  fitted  with  I 
less  than  4.  Becker's  summer  daytime  profiles  obtained  in  June/ 


July  1954. also  those  of  March/April  1958,  when  they  showed  some 


indication  of  PI,  were  fitted  with  1=4  in  order  to  see  the  dif¬ 
ference  against  1=3.  Figure  19  shows  three  examples  where  the 
fit  was  made  with  I«4  and  1=3  as  well.  Error  maps  (in  HX1  vs. 
SCI)  were  also  computed  and  are  shown  in  Figure  2o.  The  change 
between  both  I-values  is  important.  With  1=3  a  "trench"  appears 
at  rather  small  SC-values  (around  2o  km).  With  I«4,  however, 
the  same  original  profile  has  its  SC (opt)  at  much  greater  valu¬ 
es,  so  that  the  trench  is  displaced  and  is  now  almost  elongated 
along  the  SCI -axis.  Thus,  the  optimum  is  now  found  in  a  rather 
narrow  range  of  HX1 .  However,  the  error  values  show  only  little 
variation  such  that  we  have  a  large  field  in  the  map  with  al¬ 
most  equal  Err-values.  This  means  that  (when  always  determining 
the  optimum  amplitude)  there  is  a  certain  freedom  in  the  choice 
of  both,  HX1  and  SCI. 

It  should,  therefore,  be  noted  that  the  introduction  of  a  secon¬ 
dary  function  at  F-region  heights,  in  addition  to  the  main  func¬ 
tion,  considerably  changes  the  relation  between  SHI  and  SCI ,  as 
mentionned  above.  Going  from  1=3  to  1=4  makes,  in  fact,  a  con¬ 
siderable  change  in  the  reproduction  of  the  F-region.  As  the 
above  Figures  show,  a  twice- incurved  profile  can  quite  well  be 
reproduced  when  admitting  two  LAY-functions  instead  of  only  one. 
While  the  amplitude  of  the  main  function  is  always  negative,  the 
secondary  one  needs  a  positive  amplitude.  Since  the  transition 
heights  HX1 ,  HX2  were  found  to  be  not  very  much  different,  the 
effect  of  both  functions  must  partially  compensate  each  other 
in  a  certain  height  range  which  is  determined  by  the  smaller  of 
the  two  scales,  SCI,  SC2. 

Results  of  fitting  day-profiles  with  1=4  are  shown  in  Figures 
21  through  23  for  the  years  1954,  1957  and  1958,  respectively. 
The  first  of  these  years  was  at  solar  minimum,  the  last  two  were 
at  the  greatest  activity  maximum  ever  reached  in  our  lifetime. 

We  have  established  statistics  of  our  fitting  results  for  the 
two  periods  where  Becker  has  many  profiles:  summer  of  1954  and 
spring  of  1958.  The  main  results  of  this  statistics  are  shown 
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in  Table  4  which  might  be  useful  at  future  synthesis. 

Table  4.  Statistical  overview  of  the  results  of 
analysis  with  four  LAY- functions 
/Sp  =  Spring:  21.3.  -  28.4.1958; 

Su  *  Summer:  27.6.  -  1.7.1954/ 


F-region: 

HXI/km 

SCI /km 

HX1/HM 

median//quartile 

median//quartile 

median//quartile 

range 

•range  | 

range 

Sp 

272 

268-290 

9o 

— 1 

.74 

.68-. 80 

Su 

175 

152-2oo 

88 

.7o 

.62-. 82 

F-region: 

HX2/km 

j  SC2/km  | 

|  HX2/HM 

Sp 

261 

247-271 

mm 

9-21 

.69 

.68-. 73 

Su 

185 

175-209 

D 

4-48 

.75 

.72-. 81 

Valley: 

HX3/km 

SC3/km 

Sp 

117 

113-118 

7 

5-9 

Su 

114 

1o5-119 

5 

2-9 

E-region: 

HX4/km 

|  SC4/km 

Sp 

98 

9o-1o1 

5 

wm 

Su 

96 

88-1o3 

5 

WSM 

i’  teak : 


HM/km 


Sp 

Su 


371 

249 


356-392 

231-26o 


Our  results  show  that  there  is  a  clear  distinction  between  the 
main  function  (of  negative  amplitude)  which  has  an  SCI  around 
80  km,  and  the  other  one  (of  positive  amplitude)  with  an  SC2  on- 
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ly  slightly  above  1o  km.  Median  values  and  quart ile  ranges  can 
be  seen  fro®  Table  4.  When  compared  with  HM,  the  heights  HX1 , 

HX2  typically  give  a  ratio  around  o.7  instead  of  o.9,  found  in 
the  1-3  analysis: 

HX1 ,2  0.7  •  HM  (iv) 

SET  «  97  km  (v) 

With  1-3,  we  also  found  a  smaller  SCI  (63  km);  the  larger  value 
now  needed  for  SCI  (with  1-4)  must  be  attributed  to  the  partial 
compensation  of  the  two  LAY-functions  in  the  F-region. 

Like  for  1-3  we  nave  also  investigated  whether  the  SCI -values 
found  by  fitting  with  1-4  are  related  with  the  thickness  para¬ 
meters  SH  and  SQ.  No  good  correlation  could  be  found  with  SQ. 

For  SH,  however,  there  exists  some  correlation  as  can  be  seen 
from  Figure  24.  Separately  from  the  two  data  sets  we  found  the 
following  linear  relations: 

cni  fo.65  *  SH  in  spring  , 

t  o,83  •  SH  in  summer.  (vi) 

When  trying  to  combine  both  sets  (as  in  Figure  24),  we  better 
use  an  incurved  relation,  for  example  (1-4,  day): 

SCI/km  »  12o  -  0.0066  (SH/km  -  19o)2  .  (vii) 


As  shown  above,  we  found  quite  different  values  with  a  three 
function  fit.  Eoth  approaches  should  not  be  confounded  there¬ 
fore  -  except  for  the  fc-  and  V-fits  which  use  to  give  compara¬ 
ble  results.  As  for  SC2  (of  the  "correction  function"  provoking 
the  Fl-feature)  the  dispersion  of  the  points  is  larger,  the  me¬ 
dian  relation  being: 


SC2 


o.lo  •  SH 
o.o9  •  SH 


in  spring  } 
in  summer  . 


(viii) 


Thus,  the  value  of  SC2  lies  between  15  and  11sh>  of  that  of  SCI. 
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It  is  a  practical  question  whether  a  reference  reproduction 
(e.g.  in  IR1 )  needs  the  Fl-feature  to  be  reproduced.  Our  re¬ 
sults,  apparently,  show  that  one  cannot  simply  add  a  secondary 
function  in  the  F-region  description  without  admitting  some 
change  in  the  parameters  of  the  main  function.  Summarizing  our 
findings:  if  HXO,  SCO  are  the  main  layer  parameters  in  an  1=3 
description,  then  -  when  going  over  to  1*4  -  one  might  make 
new  specifications  as  follows: 

f  HX1  *  o.8o  .  HXO;  HX2  -  o.76  •  HXO 
l  SCI  *1.4  •  SCO;  SC2  *  0.12  •  SCI 

3.  RESULTS  LOWER  IONOSPHERE 

In  the  lower  ionosphere  we  encounter  a  comparable  situation  as 
we  found  in  the  middle  one.  At  the  upper  boundary  peak  height 
and  density  (now  of  the  E-layer)  are  known;  the  profile  has 
(very  rarely  a  minimum  but)  almost  always  a  turning  point  near 
8o  km.  Kechtly  and  Eilitza1  established  a  simple  empirical  rule 
for  the  day-time  electron  density  NmL  for  that  point  which  is 
used  in  IRI  in  the  following  shape: 

tNmL/m“^  =  (6.o5  ♦  o.o08  *R)  •  1o8  • exp(-o . 1 /(cos  £  )^*^) .  (15) 

h  is  the  Zuerich  sunspot  number  and  X  the  solar  zenith  angle. 

The  formula  was  obtained  with  a  combination  from  a  set  of  in- 
si  tu  rocket  measurements  obtained  with  a  combination  of  Lang¬ 
muir  probe  and  (ground-to-rocket)  phase  propagation  experiment. 
The  night  value  of  NmL)  is  fixed  in  IRI  to  4  •  1o8  m*^. 

^  E.A.  Mechtly  and  L.  Eilitza,  Models  of  D-region  Electron  Con¬ 
centration,  Rept.  IPW-WB1 ,  Fraunhofer-Institut  ftier  physika- 
lische  Weltraumforschung,  B-78oo  Freiburg  i.Br.,  F.R.G.,  1974 
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Though  the  lower  Ionosphere  was  not  the  principal  field  of  the 
present  investigation,  we  felt  that  we  should  at  least  investi¬ 
gate  whether  our  proposed  representation  with  LAY-functions  is 

2 

also  applicable  to  that  height  range  .  In  the  following  we  des- 

*5 

cribe  the  results  obtained  with  three  different  data  sets  . 

5.1  McNamara's  Data  Collection.  This  is  a  rather  complete 
collection^  of  all  available  measured  profile  data  up  to  1977. 
Most  profiles  are  incomplete  and  the  absolute  calibration  is 
often  doubtful.  Ramanamurty  has  selected  a  much  smaller  set  of 
profiles,  covering  the  whole  height  range  including  the  E-peak. 
This  latter  is  sometimes  ill  defined.  Ramanamurty  first  tried 
to  take  the  peak  density  NmE  from  empirical  formulas  describing 
the  critical  frequency  foE,  as  evaluated  on  ionograms.  However, 
the  absolute  calibration  of  the  collected  data  used  to  be  too 
bad,  so  that  the  disagreement  between  published  profile  density 
and  that  obtained  from  ionograms  was  too  large. 

Our  method  is  not  bound  to  absolute  values  -  the  description 
gives  the  logarithmic  density  ratio  to  the  peak  value.  So,  we 
work  with  relative  densities  anyway.  Therefore,  we  finally  de¬ 
cided  to  determine  an  NmE  from  the  profile  data  themselves.  This 
guess  sometimes  is  difficult,  particularly  where  more  than  one 
maximum  is  reported.  We  used  to  take  that  which  was  nearest  to 
the  standard  peak  altitude  of  IRI  (1o5  km). 


2 

An  important  part  of  this  particular  investigation  was  made  in 
cooperation  with  Y.V.  Ramanamurty  (on  leave  from  National  Phy- 
sical  Laboratory,  New-Delhi,  India).  His  stay  in  Germany  was 
sponsored  by  Deutsche  Forschungsgemeinschaft,  Bonn. 

3 

Y.V.  Ramanamurty  and  K.  Rawer,  Modelling  of  the  lower  iono¬ 
sphere  according  to  IRI  guidelines,  Adv.  Space  Res.  6  (to  ap¬ 
pear  in  1987). 

4  L.F.  McNamara,  Ionospheric  D-region  Profile  Data  Base,  A  Col¬ 
lection  of  Computer-Accessible  Experimental  Profiles  of  the  D- 
and  lower  E-region,  Rept.  UAG-67,  WDC-A  for  Solar- Terrestrial 
Physics,  NOAA,  Boulder,  Co.,  U.S.A.,  1978. 
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Analyzing  some  2o  daytime  profiles  in  terms  of  three  LAY-func- 
tions,  we  found  the  dispersion  of  the  resulting  parameters  ra¬ 
ther  large.  A  dependence  on  the  solar  zenith  angle  £  could  not 
clearly  be  established.  Since  our  aim  is  to  get  a  set  of  signi¬ 
ficant  average  parameters,  we  feel  that  the  median  values  of 
our  selection  might  be  considered  as  more  or  less  significant. 
These  are  given  in  Table  5. 

Table  5.  Median  lower  ionosphere  LAY  representation 
(for  day  conditions). 

Main  secondary  1  secondary  2 

9X  8o.4  87  7o.3  km 

SC  11  5.3  6.3  km 

A  -22.4  1.2  4.5 

This  is  probably  o.k.  for  the  first  two  lines,  i.e.  the  geome¬ 
tric  parameters.  As  for  the  amplitudes  it  might  not  be  helpful 
to  take  median  values  because  the  shape  of  the  profile  is  very 
sensible  to  the  exact  values  of  the  amplitudes.  Also,  in  the 
future  application  the  amplitudes  will  not  be  taken  from  past 
experience  but  be  determined  at  the  application  itself  (so  as 
to  satisfy  certain  constraints). 

Studying  this  in  more  detail,  we  have  computed  quite  a  number 
of  profiles  with  the  geometric  parameters  HX  and  SC,  given  above, 
but  with  different  combinations  of  amplitudes.  It  appeared  that 
rather  small  changes  of  the  amplitudes  of  the  secondary  func¬ 
tions  might  have  the  effect  that  an  otherwise  monotonuous  pro¬ 
file  is  transformed  into  one  with  a  peak  somewhere  around  or 
even  below  8o  km.  It  is  even  possible  to  "overthrow”  the  sche¬ 
dule  and  produce  a  profile  running  with  decreasing  height  to¬ 
wards  always  greater  values. 

However,  since  the  future  choice  of  amplitudes  will  be  limited 
by  strict  constraints,  it  might  be  interesting  to  show  how  this 
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choice  could  influence  the  profile.  In  Figure  25  we  show  a  few 
of  our  computed  profiles,  all  with  the  same  geometric  parame¬ 
ters.  It  can  be  seen  that,  with  an  appropriate  determination  of 
the  amplitudes,  one  can  obtain  very  reasonable  profiles. 

5.2.  Median  Profiles  from  the  Magnetic  Equator.  At  low  lati¬ 
tudes  there  exists  a  collection  of  rocket  profiles  from  Thumba 
(India),  near  the  magnetic  equator^.  All  of  these  were  obtained 
between  1966  and  1978  by  the  same  measuring  device,  the  PHL- 
Langmuir  probe.  The  profiles  cover  the  height  range  up  to  16o 
km,  starting  at  about  6o  km  by  day  and  at  about  85  km  for  twi¬ 
light  and  night.  They  are  averaged  in  four  groups  as  described 
in  Section  2. 

The  authors  give  a  mean  noon  profile  which  we  have  analyzed. 

For  the  other  groups  we  produced  one  median  profile  for  each 
group  and  analyzed  this  profile.  While  individual  profiles  show 
some  irregularities,  the  median  profiles  are  reasonably  smooth 
so  that  they  may  be  taken  as  representative.  The  E-peak  was  al¬ 
ways  matched  with  the  observed  profile  curve.  The  results  of 
this  analysis  are  shown  on  Table  6.  The  relevant  profiles  (mea¬ 
sured  and  fitted)  together  with  the  contributing  LAI-functions 
are  shown  in  Figures  26 . 


E.H.  Subbaraya,  Satya  Prakash  and  S.P.  Gupts,  Electron  Len- 
sities  in  the  Equatorial  Lower  Ionosphere  from  the  Langmuir 
Probe  Experiments  Conducted  at  Thumba  during  the  1  ears  1966- 
1978,  Scientific  Rept.  ISRO-PRL-SR-1 5-83 ,  Indian  Space  Re¬ 
search  Organisation,  Bangalore,  India,  Lee.  1983. 
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Table  6.  Low  Latitude  Irofiles  of  the  Lower  Ionosphere+ 


Type  Hours  Number  I  HX1  SCI  HX2  SC2  HX3  SC3/km 


Day 

1o4o 

12 

3 

1o1.2  13.4 

82.7  1.2 

81.8  9.5 

-H15 

(0.064) 

(0.322) 

(-4.716) 

Morning 

o54o 

1 

2 

1o1.2  6.7 

89.3  -1.7 

I 

t 

3 

(-0.293) 

(-1 . 1o2 ) 

twilight 

-o615 

1 

85.7  3.o 

(-4.170) 

Evening 

183o 

1 

1 

i  3 

3 

84.9  11.1 

92.2  8.1 

twilight 

-19oo 

! 

(-89.60) 

(38.49) 

Night 

22oo 

r — 

6 

2 

85.6  3.9 

93.5  4.6 

• 

o215 

(6.685) 

(-5.547) 

+Numbers  in  parantheses  give  the  amplitude  found  at  fitting. 


For  the  median  daytime  profile  three  LAY-functions  are  clearly 
needed.  This  profile  peaks  at  1o7  km  and  shows  a  turning  point 
around  80  km  at  a  density  of  about  1/3oo  of  NmE  /i.C.: 

log  (N/Nmax)  =  -2.5  /.  Compared  with  the  present  IRI,  this  is 
not  far  from  its  NmD/NmE  ratio.  The  lower  boundary  is  at  about 
65  km  /with  -3.5,  as  for  the  present  IRI/. 

The  twilight  flights  were  so  arranged  that  the  solar  zenith 
angle  was  between  9o  and  loo0.  At  morning  twilight  the  (here 
alone  interesting)  profile  below  HmE  can  quite  well  be  represen¬ 
ted  by  one  LAY-function  only.  As  for  evening  twilight,  the  three 
basic  profiles  are  more  different  so  that  the  median  profile 
might  not  be  representative.  So,  we  obtained  an  almost  constant 
electron  density  between  9o  km  and  the  "peak"  at  1o2  km.  At 
least,  two  LAY-functions  are  needed  ;  for  perfect  description 
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three  of  these  are  preferable  so  that  the  flat  peak  can  be  bet¬ 
ter  reproduced.  In  this  special  case  our  fit  leads  to  rather 
large  amplitudes  (because  the  particular  shape  is  brought  about 
by  compensation  of  functions  of  different  sign). 

Though  we  describe  (in  Figure  and  Table)  the  night  profile  of 
the  Indian  authors  with  the  help  of  two  LAY-functions,  this  is 
not  really  needed.  A  satisfying  description  could  be  obtained 
with  one  only. 

We  feel  that  the  results  of  this  homogenous  set  of  measurements 
should  be  given  preference  in  comparison  with  the  inhomogeneous 
set,  discussed  in  Subsection  5.1.  They  could  directly  be  used 
for  establishing  a  low  latitude  lower  ionosphere  reference  pro¬ 
file. 

It  is  hoped  that  similar  sets  of  homogeneous  measurements  might 
be  obtained  for  other  latitudes  too. 

5.3.  Median  Profiles  from  Radio  Wave  Propagation  Data. 

Singer  has  constructed  profiles  by  trial-and-error  comparison 
with  data  obtained  from  different  techniques  of  probing  the  io¬ 
nosphere  with  radio  waves.  We  analyzed  his  profiles,  mostly 
with  three  LAY-functions.  Figures  27  show  some  results. 

It  must  be  emphasized  that  Singer's  night  profile  is  rather  dif¬ 
ferent  from  the  in-situ  measured  profiles  by  the  fact  that  it 
has  a  "staircase"  aspect  with  almost  discontinuous  transitions 
from  one  level  to  the  next  one.  It  is  explained  that  a  good  re¬ 
production  of  the  radio  data  from  different  working  frequencies 
cannot  be  otherwise  obtained.  Under  these  circumstances  our  fit¬ 
ting  method  necessarily  leads  to  a  series  of  well  distinguished 
short-scale  bends,  so  that  there  is  no  difference  against  the 


W.  Singer,  Jittrchi  ...  EUKtrtfnOnkenuntratteMprofiU  •••, 

Z.S.  Meteorologie  26,  23'1-2^3  (1976) 
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original  "step-by-step"  method  of  Booker.  One  needs  more  than 
three  LAY-functions  with  this  very  particular  profile. 

We  feel  that  this  profile  should  be  considered  as  controver¬ 
sial.  The  special  features,  mentionned  above,  occur  at  very 
low  electron  density,  i.e.  at  heights  where  in-situ  measure¬ 
ments  are  very  difficult.  May-be  that  the  present-time  activi¬ 
ty  in  the  "Middle  Atmosphere  Program"  could  end-up  with  relia¬ 
ble  data  for  this  height  range  by  night.  In  the  actual  IRI, 
these  low  electron  densities  are  outside  of  the  range  of  our 
particular  interest. 

6.  CONCLUSIONS 

Summarizing  our  findings,  we  state  first  of  all  that  the  pro¬ 
posed  system  of  profile  representation  can  very  well  be  applied 
to  all  types  of  electron  density  profiles  observed  in  both, 
middle  and  lower  ionosphere. 

For  the  middle  ionosphere,  we  established  in  Chapter  4  numeri¬ 
cal  relations  depending  on  the  peak  altitude  HmF.  These  might 
provisionally  be  applied  for  mid-latitudes.  By  day,  they  give 
the  geometric  parameters  for  three  (or,  in  the  presence  of  FI, 
four)  LAY-functions .  By  night,  one  function  is  enough  for  des¬ 
cribing  the  upper  profile  part  which  refers  to  the  F-region. 

The  lower  part,  i.e.  valley  and  E-region,  must  be  inferred 
from  other  sources. 


For  the  lower  ionosphere,  in  Chapter  5  we  could  at  least  show 
that  the  description  with  three  LAY-functions  allows  a  very 
satisfying  reproduction.  From  different  data  sets  somewhat  dif¬ 
ferent  sets  of  geometric  parameters  were  obtained.  The  "Task 
Group  on  IRI"  will  have  to  decide  which  of  these  is  to  be  pre¬ 
ferred  -  or  whether  new  basic  data  should  be  looked  for.  At 
night,  the  lowest  electron  density  features  at  the  base  of  the 
ionosphere  (Singer's  night  profile)  might  provisionally  be  ais- 


regarded.  Then,  one  LAY-function  describing  the  Hnite-E" 
bottomside  is  good  enough. 

Finally,  future  application  of  our  proposed  system  in  Ihl 
necessitates  the  establishment  of  suitable  geophysical  con¬ 
straints  from  which  the  amplitudes  should  be  determined  indi¬ 
vidually,  i.e.  in  the  IRI  computer  program  itself.  These 
should  not  arbitrarily  be  distributed  over  the  height  range 
but  must  be  chosen  in  a  way,  that  to  each  LAY-function  at  least 
one  relevant  condition  exists. 


CAPTIONS  TO  THE  FIGURES 


Fig.  1.  LAY-function  of  altitude  z  (ordinate),  variation  of 
geometric  parameters: 

(a)  SC  fixed  to  1o  km,  HX  varied; 

(b)  HX  fixed  to  25o  km,  SC  varied. 

Fig.  2.  Error  map,  i.e.  reduced  error  sum  as  function  of  HX 
(ordinate)  and  SC  (abscissa).  Arrows  show  directions 
of  "steepest  descent". 

Fig.  3.  LAY-function  multiplied  with  amplitude  A.  Interdepen¬ 
dence  of  A  and  HX  for  fixed  SC.  Amplitudes  were  so 
chosen  that  same  curvature  occurs  at  peak  HM. 

Fig.  4.  Flux  diagrams  of  optimization  procedure  APTR: 

(a)  Main  program;  (b)  1-6  Subroutines. 

Fig.  5.  Printouts  of  final  optimization  program  APTR  (in  SIMULA 
computer  language): 

(a)  Main  program  APTR; 

(■^^Subroutines  used  in  APTR; 

(L3)  Subroutine  LSFURR  followed  by  CHOLESKY  (solving 
a  system  of  linear  equations). 

Fig.  6.  Night  profiles:  (a)  normal;  (b)  odd. 

Fig.  7.  Error  maps  (see  Figure  2)  to  Figure  6 . 

Fig.  8.  Error  maps  for  eight  different  night  profiles  of  1954, 
1957  and  1958.  (date  and  hour  see  lower  left-hand 
corner). 

Fig.  9.  Ratio  (as  function  of  HK)  of  greatest  admissible  HX 
(from  night-time  error  map),  HX1 ,  to  (a)  HK,  (b)  HHA, 

(c)  Hl<U.  Thick  line:  connection  between  median  points 
of  1954  and  1958.  Thin  lines:  delimit  quartile  ranges. 


t  t  , 
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Fig.  1o.  Difference  (as  function  of  HM)  between  greatest  ad¬ 
missible  HX  (from  error  map),  HX1 ,  and  (a)  HM, 

(b)  HHA,  (c)  HvJU.  /See  caption  to  Figure  9/. 

Fig.  11.  Ratio  of  half  density  height,  HHA,  and  peak  height, 
HM,  as  function  of  HM.  /See  caption  to  Figure  9; 
o  1954,  x  1957,  +  1958/. 

Fig.  12.  Night-time  plasma  frequency  profiles  (right  hand 

side)  and  cuts  through  the  error  map.  (one  LAY-func- 
tion  only)  at  HX1  *  HM/2  (full  line)  and  at 
HX1  -  HM/3  (broken  line). 

/Becker1 s  profiles  from  00  to  04h  on  25  June  1954; 
hour  indicated  inside  profile  drawing/. 

Fig.  13.  Cuts  through  error  maps  (one  LAY-function)  at 

HX1  -  HM/2  for  all  night  profiles  given  by  Becker  for 
26/27  June  1954.  /Hour  in  the  left-  or  right-hand  up¬ 
per  corner/.  In  Subfigures  (a),(b),(c)  the  cuts  are 
arranged  after  the  value  of  SC  at  minimum. 

Fig.  14.  Survey  to  Figures  13  showing  the  variations  of  the 
optimum  (i.e.  position  *  SC  for  minimum, and  minimum 
Err-value,  Err  min).  Further,  the  width  of  the  mini¬ 
mum:  Del ta^ between  2*Err  min  points,  and  delta2  at 
Err  »  5o.  HX1  «  HM/2  in  bottom  line/(night  26/27 
June  1954,  18-06h)/. 

Fig.  15.  Correlograms  between  optimum  SC,  abscissa,  and  lower 
half-density  layer  thickness,  SH  (top),  or  lower 
quarter-density  thickness,  SQ  (bottom).  Night  profi¬ 
les  from  1954  (solar  activity  minimum)  to  1958  (maxi¬ 
mum). 


■O.* 
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Pig.  16.  Reproducing  day  profiles1  with  three  LAY-functions : 

(a)  Reference  profile  with  an  E-P  valley; 

(b)  Eecker’s  profile  886  with  no  FI  indicated; 

(c)  Becker's  profile  74  with  PI  indicated  (but 
overlooked  in  the  fitting  procedure). 


Fig.  17.  Six  day  profiles  without  FI  fitted  with  three  LA1- 
functions  (1=3).  Pecker  profiles  numbers:  541,  661, 
761,  826,  953,  1325. 

/G  »  Gulyaeva's  estimate  of  HHA/. 

Pig.  18.  Error  maps:  isolines  of  reduced  error  sum  in  HX 

vs.  SC  coordinates  of  main  LAY- functions.  Example 
for  day-time  profile  without  FI  (1=3). 

(a)  Iiarge  scale  survey; 

(b)  Neighbourhood  of  minimum,  small  scale. 

Fig.  19.  Three  profiles  with  pi-features  fitted  with  1=3 
(a,c,e)  and  1=4  (b,d,f)  for  comparison.  Becker 
profiles  numbers  74  (a,b),  114  (c,d),  1484  (e,f). 


Pig.  2o.  Error  maps  for  three  profiles  with  Fl-feature, 
fitted  with  1=3  (a,c,e)  and  with  1=4  (b,d,f). 
Same  Becker  profiles  as  in  Figure  19. 


In  our  profile  drawings  the  original  data  points  are  drawn 
as  circles.  The  individual  contributions  of  the  different 
LAY-functions,  A  •  LAY,  are  shown  with  different  symbols. 
Their  sum,  i.e.  the  approximation  finally  reached,  is  given 
by  triangles A.Note  that  the  drawing  program  gives  straight 
line  connections  such  that  the  marked  points  only  should  be 
considered. 
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Pig.  21. 2  Six  day-time  profiles  with  Pi-feature  fitted  with 
four  LAY- functions  (1-4):  Becker  profiles  numbers: 
122,  177,  18o,  185  (two  versions),  294,  3o7. 

All  of  1954  (solar  minimum). 

Fig.  22.  Bay-time  profile  of  1957  (Becker  number  541,  solar 
maximum)  fitted  with  four  LAY-functions  (1*4). 

Fig.  23.  15  day-time  profiles  of  1958  (solar  maximum)  with 

four  LAY-functions  (1*4).  Becker  profiles  numbers: 
661,  761,  821,  826,  953,  957,  11o4,  1223,  1275 *  1325, 
1365,  1428,  1545,  1613,  166o. 

Fig.  24.  Correlograms  between  the  optimum  scale  SC1  of  the 
main  LAY-function  (ordinate)  and  the  lower  half¬ 
density  thickness  SH  (abscissa),  obtained  from 
day-time  Becker-profiles  fitted  with  1*4. 

Fig.  25.  Lower  ionosphere  profiles  (day)  computed  with  three 
LAY-functions  (1*4)  using  the  median  geometric  para¬ 
meters  obtained  by  analysis  of  McNamara  profiles 
(Table  5)  and  different  amplitude  combinations. 

(a) -A-j  *  parameter,  A2  *  5.2,  Aj  *  4.3; 

(b)  A1  *  -22.4,  A2  *  parameter,  A^  -  4.5; 

(c)  A1  *  -22.4,  A2  *  4.8,  A?  =  parameter. 


2  Since  the  representative  function  is  fitted  to  the  input 
points  (circles),  it  might  show  some  deviations  in  between, 
which  are  not  necessarily  shown  by  the  measured  profile. 

The  two  versions  of  profile  185  show  what  might  happen  in 
the  E-F  valley  when  the  height  difference  of  the  data  points 
becomes  too  large.  This  can  be  avoided  by  introducing  inter¬ 
mediate  data  input  at  these  critical  intervals. 


»■  1  k.  <rVir .X.  iA  A  WWL.V-*  AVLV.'Jl?, 
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Fig.  26.  Median  profiles  after  rocket  measurements  at  Thumba 
(India  -  magnetic  equator,  see  Table  6).  Standard 
analysis  with  up  to  three  LAY-functions  for: 

No.  1  »  noon;  2  »  night;  3  *  morning  twilight. 

No.  5  is  the  same  night  profile  as  No.  2  analyzed, 
however,  with  the  upper  peak  appearing  around  16o  km. 
/This  example  shows  that  our  analysis  is  also  suited 
for  abnormal  conditions/. 

Fig.  27.  Analysis  of  standard  ionosphere  profiles  as  esta¬ 
blished  by  W .  Singer. 

(1-4)  Summer,  Covington  index  F  =  75;  solar  zenith 
angle  (sza)J(l)  8o°,  (3)  6o°,  (4)  4o°. 

(5-8)  Winter,  F  «  75;  sza:(5)  85°,  (6)  8o°  average 
quiet  conditions,  (8)  8o°  low  absorption. 

(1o)  Summer,  F  «  1 5o ;  sza  7o°. 

(13-15)  April,  F  =  75;  sza;(13)  8o°,  (15)  6o°. 

(16)  F  =  75;  -  night  conditions  -. 
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f  MCI*  IIUUI  l| 

J  UlHfiH; 

4  ft*  l»M  lfl>  1  Wftl  «M  M  MtUlt*  If  Ukl(tl)  ta  t  l*|» 
f  fli)«Uili)<lk  (U(  f(i)f  (ffulit  at  a  iiht  a«T|»; 

f  fH  H't.i  »•  U*f  ftlKI 

t  IM  i  (IMHIt  I  Ml  I  »  I  f  T  I  , 

•  MMCW1I  lailTlj  (IMtkT  UIM  l|T«M  it  Mill  Mf, 

98  Mala  utMl  u; 
ff  aCs*ti«f{ 

fa  im  at*  1,2  »•  kmiiiwmi'Dj 

fa  IM  lt>f  fill  t  f«IU  M4  M  MitUll  l«  I  fat* 

fa  ia>«!*Mim*N  it.sc  Mil,  tfanit  at*  staff  aftt#* 

fS  UltM(M|flM|l^il.l,f,llll; 

U  mflf»tM|MMIII.f|l,ltf,M^fCiM.,ii;i;  till  Mil* 

It  IMi  COXtaf  IMS  I  a  t  f  T  t  i 
fa 

f*  fHUMM  IMItlM  ^iM*Mla (Mil  ,aa,M (l,t ,a .till  ,*l( I ; 
fa  (mhii  imi  im  tt*aiMCf  to*  M  iiuntt  ttctafiat  aa-atciaa  caia,aa 
ft  taaa.aa  •>  m*i«mmi  im  faacitM  t««  m  an  imhuii  •*  Mtu*iiM 
it  Mk  at  •  aciaaaaaa  Mini  cata  taiifti  a*  mmimi  .hmi  m  taaaiMitt. 
fa  if  aa  ifMiuaui  vaaiafiaa  ta  mat  tii  j«m  ta  u*u  *ai. 
fa  UtiCMM  a#  m  mufiii  ant  a*  unit  aaaa«  x(ft*»>  at  itrttittaatt 
ff  Mai  Coot  a#  aaat.  (Mtu>,  aaaaat  as  ia#»ft  ff  m . 

fa  M«tau,  aa  latimui  »•«  imi, 

ft  aaat  ctiij  ititMl  iim*.  mm  xa,  lavcaa*  taaat  taaf ; 

fa  aaaat  t»ia,aa.»,«a.a.T,  aaaa  MMImM  i#a.  a*a«t  atn 

ff  aaata  lartfaa  m.m^i,  mm  ii.mmj 

If  mkiM,  tOMWaf  MU. 

ft  aaata  aaaav  ta*McCi4»),a»<f taai.aaif ,xi. 

ff  faa  tin  sr«a  l  ntu  at  0#  MUHHadi, 

9S  If (tl  i*lf IIMMlMIMtSMif tal . 

fa  aiMt  m  n»i  staa  «  titii  aa  m 
ff  if  saac at  ta  t  mii  aaata 
fa  k»M*1,  MlllitlMKHMIi 

ft  MMliMIKIItlltMIMItlflMl, 

fa  aataioaaui-aniaf  faa  tatt  af<a»i>*itai,  aat«a. 

f«  atiHdM  4i.»  IMI  I  k*«lk  M  M  I'  lulUl  la  I  IX* 

*4  atata  a»*t*a»iMM<ai>«i<aiifa<i«atj 
at  »i««4,itt  t*4  i 

42  aa>«aaafi«al, 

41  If  aa  i.T  4-f  'Ml  MU  IK, 

44  IM  ||«1  ttta  1  aatn  M  M  |l  IIUUI  It  I  fat* 

41  tta<tli(*l !•**«« ttmm,  tat. 

4a  laaatat  uMi  i  taac  •  t  i  l  a  «aa, 

41 


4«  tatt  aaactaaac  •iaaaitj,«a.ca,ftiai . 

4t  i aa*ta i  ai*«t  a*  •him*  «  m*Mi*  'mim*  l  miimmi 


ts  ti  fantaittaai  mkuui  it 

If  t*  rkHiiMi  lataaxat  I»intti,i4t4t 
$•  tax  nmnum  M,*lh  .a.aaatt  II, 
tl  taata  fait. 


tax  tain. 
HI.  tl. 


ft  Mat*  tatt  Cl.lt, 

•t  aat>|f  till!  t*  »*n»  T««a  <mi 
fa  tti*a.t*t*«ti-tTi;i44.t4,i4fi4 
t»  n  i«  a*  ta  ta«a  ta*t  »*n. 
ta 

•a  aiattf-i ’/«*•<?!,  *•*. 

M  c  aaata  i  taaf  a  * 


till  fan, 


# 


* 


•a 
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5  b* 


(2  IIU  MtCIIUll  kAY<41,aC,ltl6),  (674(67  0«€  *07,70047. 4. , SC, 

4)  J.UJT  A69kl1U6t,WITN  kC((..A1>*3  Al  I*.  «»«»  Juk.  *1, 

4*  IUI  NI.S(«l|Hj  M«» 

41  4,*f i moyou, sc.ud-c 7Ta(*9,sc,„i i -<«-**»•  (7Stc*»,sc.h* i /$c  roo; 

44  (MMIt  t«H  LAY. 

44  MAI  700C(6UAf  AY  UIA  Clli  ,4k  ,»9  A  ,  7(  ,  M  ,  Y  ..  )  ; 

44  (MfiYT  IA.fl  YUM  tS.*70.iU6,0((BuC16,9(  kllUOO  C*0  >(>007  700  077I4U7); 

76  I4TIMA  UL. Ml;  ACAL  MA; 

71  AOOAf  44...I.Y; 

72  um  14T((fO  ((,«(;  AfAL  7i.Y7.iu; 

71  AA4AV  079(1:41) ,7<1:9k) 

7*  hi t *2*41;  ^ 

74  17  90A7AU1  7414  61614 

77  6«ftM*f;  OuTT(*T(-  (  Y->; 

76  76.  *0t*1  17(7  1  U4T1L  "k  BO  0(814 

7*  tuYlMII,  tl*U*l(ll);  047147(4.10); 

04  ■t*4(*V(«A>;  0UT14Y (I .16)  140,  (40; 

•1  6Ul*7!i*6;  760  c«:*1  11(7  1  4.471*  41  BO 

•2  *0414  I7i*f(u)*f(l«l;  7l«.»i**»*Y»*4(*«);  1u:*Su**Uk>  (40; 

91  A7uaoi*Yi/sw; 

6*  (644661  (461  A  7  U  0  0;  (40; 

•1 

64  A04L  700( (BW0(  (774(1.41.11). 

07  (647(47  AA.(A  *A»  02.(711(14  7(7*4017104  1A04  0  A7  L(7Y 
•6  10  l(6(A(  .174  1(071  1/01  Al  B1(4T.  .1  *WlT  •(  46  0; 

60  40*1  >.6l. «t. 

46  oi.i4  ((At  o1;oli*(i>il)/61; 

«1  17  *6(1011  61  .0  74(4  4(614 

92  (710  i*l  7  al  67  C  74(4  .1  lltl  .  (40  HU 
01  (770 i *(44 1 *((7101) )  146; 

•4  (044(47  ( 96f  (77  'i 

90  4046  746C(0VA(  (71 7 1 1 .61 . « 1 1 ; 

94  (044(41  04. (A  44Y  02. 

07  (9.7(16  1**7  7446  C  AY  11*7  10  1  07  01647  .  (711 *01 01*01 1 ¥ (  07  (77*. 

06  41  76(7  Of  4(  A; 

99  6041  A. 61.01, 

100  600(4  ((At.  al;  .1  i*(l-01>/01; 

1«1  17  60(191)  67  AC  74(4 

lot  6(614  (7(11*17  ('  (1  C  74(4  1  Ik  16  3  (49  fkfi 

16?  (917 |*1|(1*6A6C«611)  440; 

164  (694047  (641  (  7  1 

16? 

1(6  Off Ak  766(101741  (MAltld. >1.90.760, Al, *01.1. I... 1I0>; 

1(7  (6*4(67  l A *0 A  *AV /i.l.*A.II *0| B.( t B  1(0(0 |U* (  1  ..* 1 ) **»I0A6(  10WAA6B  Bl«. 
106  41  46.07  .  A  V  *  I  (  A>  1  .174  464*2*41  *(0*  .(0(7  7  .761 1  t  *(6 )  .  7  *1U4(A49L  11*L  07  )• 
100  149.1.  0(0411  4A440I66  697 I *.9  (444171(4  *1  7*6, *3), 110. 

116  9(16>1  0(64(14  4466(0  (OW6t, 

111  44* t  *61.734. HA;  167*0*0  *1,41.  A«*k  11?.  *00*1  A 1 ,*0? ,7 .74 ,704 , 4. 7; 
lit  60(14  1671MB  11. **,444.  0(A.  ft, 11. 1U, 77;  6(41*1*41; 

11*  66  i  *7  i  :  *C ;  7(0  1  11(7  1  .471.  >1  B«  •((!•>  1 I : *7  I «•  > -7  <*«  ) ; 

114  *1 I*71*«l*4)*t7*t 7;  tk:*1«**«*B>  *64; 

111  4(0411 r*ll t*7| 7*y, 


114  17  *1  It  lie  74*6 

117  7.0  OOt**  1117  1  66 7 1 k  4 7  14  A C .»*•):*• 1  lot ) ; 
lit  '**  1 1  :  •  1  till  1  y411i  6(6  B(  9.4111)1*74(111 
119  ( **M  4 1  (64*  4(0 


4(416  110 1*11 ; 


*M, 


*. 


(60; 
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5  b-3 


1  PROCEDURE  LSFURRtM.PL,MPA,PL,AHPL,X,Y,.,F>; 

2  COMMENT  RAWER  JUNE  65  WITH  GIVEN  COEFF  AAR  A  T  PLd...2*M>  AND  N1 

3  AiAttkS  IN  REPRESENTATIVE  FUNCTION  F<1..*l)  FINDS  BT  LEASTSSUARES 

A  (CNCLESkT)  OPT  I  HUN  AMPLITUDES  AMPu(1...Nl»  SO  THAT  . El 6HT ED t 1 . "Li 

9  RREDUCED  SQUARED  ERROR  SUN  COF  F  AGAINST  MEASURED  VALUES  V?  IS  MINI* 
3  MIXED.  THEN  FUNCTION  tARRAT  Fl  COPPUTEO  WITH  THESE  AMPLITUDES. 

r  PROCEDURE  CHOLESKY  t  REAL  PROCEDURE  LAY  NEEDED; 
t  NAME  AHPLtF;  INTEGER  N1.Nl;  REAl  HNA; 

*  ARRAY  AMPL iF.PL.W.X.Y; 

10  BEGIN  INTEGER  1.J.N.L.N6.  REAL  F IL • FU AI J , FU IL • P2l . P2IH; 

11  N6:S2*NI;  BEGIN  CONPENT  BLOCK  ; 

12  ARRAY  FF1(1:KI.1:ML)«FUA<1:NI.1:N1),FUB(1:N1>,PU(1 SNG) ; 

13  FOR  J  :«1  STEP  1  UNTIL  N6  DO  PU < J > !*PL Cj > ; 

14  FOR  11-1  STEP  1  UNTIL  Hi  DO  BEGIN 

15  Ks*2«I;P2I:»PU(K>;P2IMs»PuU-1>; 


1*  FUEL  :«0;  FOR  L:«1  STEP  1  UNTIL  Nl  DO  BEGIN 
17  FIL:*FFI (I ,l) :«LAY(P2lP,P2I,*tL>.H"A); 

IB  FUIL :«FUIL*Y<L)*F1L«NCL) 

19  FUG(1):>FU1L 

20  FOR  1 :«1  STEP  1  UNTIL  Nl  DO 

21  FOR  JS«1  STEP  1  UNTIL  Hi  DO 

22  FUA1 J :«3; 

23  FOR  l:«1  STEP  1  UNTIL  NL  BO 

24  FUAIJ  :*FUA1J«FFI(1 ,L ) *F F I ( J ,u ) «NiL > 

25  FUA(J,1):>FLA(1,J) :>FUAXw 

24 

27  CMOLESRYd ,Nl .FUA.FUG.A-PL.SlNfaA); 

20  GOTO  CHOlA. 

29  SInGA:  OUTTE *T (“  ANPl.SIhG. 

3u  CHOLA:  for  LS*1  step  1  UNTIL  NL  DO 

31  f<l);«0;  for  j:-1  step  1  until  ni  do 

32  FCl)  :"F(l)*ahPl(j)*FFI(J,L) 

33  CONPENT  BLOCK  t  ENDE  L  S  F  U 


END; 

END; 

BE6IN 

BEGIN 

BEGIN 

end; 

END; 

END; 


BEGIN 

END;  END; 

R  R;  END; 


34 

33  PROCEDURE  CHOLESKY  tnl,  n2.  A.  B.  «.  SING),  COMMENT  RESOLVES  LINEAR 

34  SYSTEN  OF  EOS.,  EXCEPT  .HEN  SINGUlARY, 

37  NAPE  I;  INTEGER  Nl,  N2,  ARRAY  A,  4,  Y;  LABEL  SING; 

34  BEGIN  INTEGER  1C.JC.kC,  REAL  H; 

39  FOR  !C:-N1  STEP  1  UNTIL  n2  DO 

Ail  BEGIN  X(IC):*B(IC);  FOR  ;C:«1C-1  STEP  -1  UNTIL  Nl  DO 

41  BEGIN  h:«A«JC,JC )*A(1C ,JC). 

42  A(IC.IC) :»A(IC»I!)“H*A(iC.JC),  I(|Cl:*»(IC)*A(lCiJ{)*l(JC); 

43  FOR  KC:>IC*1  STEP  1  UnTil  <U  GO  A (KC . IC )  : «A (K C , I C ) -n» A < KC , JC ) 

44  END;  IF  a(IC.Il)  LE  O.u  ThEn  «0T0  SING; 

45  FOR  KC  : • I C ♦ 1  STEP  1  UNTIL  *<.  vO  A(KC.IC)  :»AIKC,1C)/A(IC,IC>  END  ; 

44  FOR  !C:«Nc  STEP  -1  UNTIL  S’  DO  BEGIN  X  ( I C >  :  »I  1 1 C > /» < I C . 1 C ) ; 

47  FOR  KCS*IC*1  STc P  1  UNTIL  Ni  DO  X ( I C ) !«X (I C ) -A <«C , I C ) •* U C )  END; 

43  END;  COMMENT  ENDE  ChuLESkY; 


yy.yy.v.y.yy 


■1.0 

L0G(N/Ntim() 


•  * 


.  net GMT 

7»  »j  1*3  Mjl'JItOITCIr',  lAIjjOJIOijO./Tuitf 
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BECKER  PROFILE  1223 


23 


ism 


ttCISMT  i  I  ihCttttT 

l  u  110  |l»  130  IM  IU  SO  ts  70  7*  M  |  Vs  *> 


2 


SUBBfiRfltfl  PROP  [  LC  3 


•13  -1.1 

LOC(N/fMtX) 


